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ABSTRACT 

Pronounced  spatial  acoustic  interference  patterns  versus  horizon¬ 
tal  range  were  observed  near  the  deep  ocean  bottom  from  a  220-Hz  CW 
point  source  also  located  near  the  bottom.  This  Is  a  result  of  the 
direct  water  wave  field  interacting  with  a  strongly  reflected  field  from 
the  ocean  sediments.  If  the  bottom  structures  are  horizontally  strat¬ 
ified  media,  the  reflected  field  can  be  expressed  by  the  Sommer feld 
integral  whose  Integrand  is  weighted  by  a  total  plane  wave  reflection 
coefficient. 

In  this  thesis  the  reflected  fields  are  analytically  evaluated  in 
order  to  understand  the  physical  mechanism  of  the  propagating  waves 
through  the  bottom  layer  structure.  Based  on  a  canonical  geoacoustic 
model  ,  the  Langer  uniform  asymptotic  approximation  is  applied  to  evalu¬ 
ate  the  total  plane  wave  reflection  coefficient  accurately.  Then  the 
Sommerfeld  integral  is  evaluated  by  the  method  of  steepest  descents. 

The  expression  for  the  reflected  field  is  fully  analytic  and  provides 
physical  Insight  into  the  understanding  of  energy  partitioning  in  the 
sediments.  With  inclusion  of  the  new  method  to  analytically  express  a 
sound  velocity  profile  shape  based  on  measurements,  the  approximate 
solution  provides  a  full  analytic  expression  of  the  reflected  field 
based  on  a  realistic  geoacoustic  model. 

Various  geoacoustic  models,  with  emphasis  on  changes  in  the  shape 
of  the  sound  velocity  profile  in  the  sediment  layer,  are  examined  and 
their  parametric  dependences  on  the  associated  reflected  fields  are  in¬ 
vestigated.  The  conditions  neccessary  for  the  shape  of  the  profile  to 
be  Important  are  also  quantitatively  defined  and  are  demonstrated  in 
theoretical  data  as  well  as  in  real  data. 
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Chapter  1:  Introduction 


tt  has  been  observed  In  the  deep  ocean  that  the  total  acoustic 
field  generated  by  a  harmonic  point  source  placed  In  the  water  column 
near  the  ocean  bottom  Is  characterized  by  strong  Interference  between 
the  field  reflected  off  the  bottom  and  the  direct  field  throucji  the 
water  column.  The  amplitude  of  this  total  field  thus  shows  a  distinct 
spatial  interference  pattern  with  horizontal  range (cf.  Fig.  1.1). 

This  thesis  involves  the  analytic  computation  of  the  bottom  re¬ 
flected  acoustic  field  using  assumed  geoacoustic  models  in  order  to 
understand  the  mechanism  of  propagating  waves  in  ocean  bottom  sedi¬ 
ments.  The  analytic  results  obtained  here  are  based  on  asymptotic  ap¬ 
proximations  valid  In  the  far  field  and  characterize  different  types  of 
rays  which  contribute  to  the  total  reflected  field.  Therefore,  they 
provide  physical  insight  into  the  understanding  of  energy  partitioning 
in  the  sediments. 

The  methods  developed  in  this  thesis  will  be  applied  both  to  the 
examination  of  the  fields  generated  by  various  "canonical"  geoacoustic 
models  and  to  actual  experimental  pressure  field  measurements  made  by 
Frisk  and  his  co-workers [5].  Extensive  tests  to  determine  whi ch  major 
ray  paths  contribute  to  the  total  reflected  field  will  be  performed 
throughout  horizontal  range  as  well  as  an  investigation  of  the  effects 
of  geoacoustic  parameters  on  such  rays. 


Fig.  1.1  A  typical  pressure  field  measurement  (Ref.  5). 


1.1  Background 


1.1.1  Oblique  Incidence  bottom  acoustic  experiments 

A  typical  experimental  configuration  by  Frisk  et.  al  .[5]  is  shown 
in  Fig.  1.2.  A  220-Hz  CW  source  is  suspended  on  the  trawl  wire  of  the 
ship  vhich  slowly  arl  fts  away  from  a  receiving  system  that  consists  of 
two  receivers  moored  on  the  bottom.  The  source  is  maintained  at  about 
100-200  m  above  the  bottom  as  it  moves  away  from  the  receivers  out  to  a 
horizontal  range  of  a  few  kilometers.  By  sampling  the  incoming  signals 
at  a  sufficient  rate.  Frisk  obtains  spatially  continuous  acoustic  fields 
near  the  bottom  due  to  the  time  harmonic  point  source  which  is  also 
placed  near  the  bottom. 

Since  the  source  and  receiver  are  placed  relatively  close  to  the 
bottom  in  the  experiments  as  shown  in  Fig.  1.2,  except  at  relatively 
short  ranges,  most  of  the  rays  interact  at  small  grazing  angles  with  the 
sediment  layer  with  increasing  sound  velocity(cf.  Fig.  1.3).  This  is 
not  the  case  if  the  source  and/or  receiver  are  placed  at  a  high  eleva¬ 
tion  above  the  bottom,  for  example  if  both  are  at  the  surface  or  if  the 
source  is  at  the  surface  and  the  receiver  at  the  bottom  as  done  in  more 
conventional  deep-sea  seismic  experiments  for  ranges  up  to  a  few  kilo- 
meters(cf .  Fig.  1 .4) . 

Spudich  and  Hel liter ger[ll ]  computed  the  seismic  responses  of  simi¬ 
lar  ocean  bottom  models  covered  by  a  4-km  thick  water  column  in  a  fre¬ 
quency  band  of  1-30  Hz  for  both  source  and  receiver  located  near  the 
ocean  surface.  Step-function  source  responses  were  computed  over  the 
range  of  2-20  km  with  the  increment  of  1  km.  In  their  computation  they 


Fig.  1.2  A  schematic  of  the  experimental  configuration  {Ref.  5) 


applied  the  generalized  ray  theoryCl]  to  major  ray  paths  excluding  any 
Internal  multiple  reflections  and  absorption  In  the  sediment  layers,  and 
concluded  that  the  sediment  layer  gradient  and  Its  thickness  were  Insig¬ 
nificant  to  the  seismic  responses.  A  major  source  of  the  effects  they 
studied  was  the  basement  structure  under  the  sediment  layers.  In  these 
cases  the  rays  interact  with  the  sediments  at  steep  angles  and  are 
therefore  relatively  Insensitive  to  the  details  of  the  sediment  velocity 
profile. 

Because  of  these  contrasting  bottom  acoustic  Interactions  due  to 
the  different  relative  locations  of  the  source  and  receiver  with  respect 
to  the  bottom  as  shown  in  Figs.  1.3  and  1.4,  we  shall  not  be  concerned 
with  geoacoustic  models  with  multi-layered  structures,  but  Instead  shall 
concentrate  on  the  detailed  profile  shapes  of  a  single  layer.  Also,  for 
higher  frequencies  used  In  Frisk's  experiments,  the  effect  of  absorption 
In  the  sediment  layers  now  become  more  important,  and  the  number  of 
major  ray  paths  may  be  greater  for  the  close  range  fields  where,  for  ex¬ 
ample,  the  effects  of  multiples  may  not  be  negligible. 

V1dmar[13]  and  R.  Chapman[3]  have  reported  that  for  typical  ocean 
bottom  structures  there  is  a  minor  energy-loss  mechanism  due  to  the  pre¬ 
sence  of  shear  waves  in  the  bottom  for  near  grazing  incidence.  Also,  In 
our  experience  with  real  data  analysis,  there  has  been  success  with  geo¬ 
acoustic  models  that  exclude  shear  waves[5],  and  we  shall  therefore  not 
Include  shear  In  our  geoacoustic  models. 


1.1.1  Sommerfeld  Integral  and  geoacoustic  models  of  the  bottom 
For  harmonic  time  dependence (e'*ut),  the  reflected  spatial  part 
of  the  acoustic  pressure  field  P  In  a  homogeneous  medium  due  to  a  point 
source  placed  In  that  medium,  which  lies  over  a  second  horizontally 
stratified  medium  can  be  expressed  by  the  Somnerfeld  IntegralCl  ,2] 


where  p  *  horizontal  slowness  (ray  parameter) 

Zq,  z  -  source  and  receiver  heights  above  the  boundary 
r  ■  horizontal  range 

R(p ,u)  *  plane  wave  reflection  coefficient 
yq  a  vertical  slowness  in  the  source -receiver  medium 
Hq^  a  the  zero-th  order  Hankel  function  of  the  first  kind 
u  «  angular  frequency  of  the  source 

Given  the  small  sound  velocity  gradient  In  the  water  column  In  the 
deep  ocean  (typically  less  than  .016  s_1ClO]),  It  is  reasonable  to 
assume  that  the  water  column  Is  homogeneous,  especially  for  experiments 
which  involve  small  source  and  receiver  heights.  In  many  abyssal  plain 
regions,  one  may  also  assume  that  the  ocean  bottom  is  horizontally 
stratified.  In  such  cases,  Eq.  (1.1)  Is  a  valid  expression  for  the 
bottom  reflected  acoustic  fields. 

The  plane  wave  reflection  coefficient  R,  which  is  the  weicjiting 
factor  In  the  Integrand  of  the  Sommerfeld  integral,  contains  all  of  the 


acoustic  properties  of  the  bottom  as  well  as  satisfying  the  boundary 
condition  at  the  ocean  bottom  Interface.  The  physical  and  mathematical 
complexity  of  R  depends  on  the  detailed  geoacoustic  properties  of  the 
bottom  structure. 

In  order  to  obtain  analytic  solutions  for  P  in  Eq.  (1.1),  we  must 
consider  reasonably  simple  geoacoustic  models.  Yet  such  models  must  be 
realistic  enough  to  generate  fields  which  compare  favorably  to  real 
oata.  hacpherson  and  Frisk[5]  have  obtained  analytic  expressions  for 
the  total  field  using  the  method  of  normal  modes,  but  their  geoacoustic 
model  was  a  slab  configuration,  a  homogeneous  layer  between  two  homo¬ 
geneous  halfspaces.  Using  a  modified  uniform  asymptotic  method. 

Stickler  and  Ahluwal  1a[12]  obtained  analytic  solutions  based  on  a  more 
realistic  model  containing  a  layer  wl  th  a  l/c2-l  inear  profile  between 
two  homogeneous  halfspaces.  However,  the  sound  velocities  were  required 
to  be  continuous  at  the  layer  boundaries,  i.e.  no  discontinuities  were 
allowed.  In  both  cases,  the  geoacoustic  models  used  as  input  to  them 
lack  some  realistic  features  necessary  to  discuss  real  data.  A  more 
complete  geoacoustic  model  should  contain  an  Inhomogeneous  layer  with  a 
monotoni cally  Increasing  sound  velocity  profile  bounded  by  homogeneous 
halfspaces  with  arbitrary  boundary  conditions  as  shown  in  Fig.  1.5. 

1.2  Asymptotic  approximation  of  the  Sommerfeld  Integral 

Our  goal  is  to  evaluate  the  Sommerfeld  integral  analytically.  This 
Involves  two  steps.  The  first  step  Is  to  obtain  the  reflection  coeffi¬ 
cient  R  based  on  our  geoacoustic  models  and  the  second  is  to  evaluate 
the  Integral.  In  order  to  obtain  R,  we  must  evaluate  the  field  inside 


A  geoacoustic  model  considered  in  this  thesis.  The  sound 
velocity  In  the  layer  monotonically  increases  with  depth. 


the  layer.  Provided  that  the  sound  velocity  within  the  layer  monotonic- 
ally  Increases  with  depth,  we  can  apply  the  Langer  uniform  asymptotic 
approximation( L.U.A.A.  )[7 ,  8]  to  compute  the  field.  The  method  has  been 
applied  by  Cormier  and  Richards[4],  and  Kennett  and  Illingworth[6]  in 
numerical  seismic  synthesis.  The  advantage  of  the  L.U.A.A.  is  that, 
unlike  the  W.K.B.J.  method,  It  is  uniformly  valid  for  any  horizontal 
slowness  (or  equivalently  horizontal  component  of  the  wave  nunber)  in  a 
given  medium  for  any  profile  shape  as  long  as  the  sound  velocity  mono- 
tonlcally  increases.  Once  the  field  Inside  the  layer  is  obtained  analy¬ 
tically,  the  analytic  expression  for  the  reflection  coefficient  R  can  be 
founa  easily  from  the  boundary  conditions. 

It  is  of  interest  to  understand  the  energy  partitioning  of  pro¬ 
pagating  acoustic  waves  associated  with  various  ray  types  as  a  function 
of  horizontal  range.  In  order  to  do  this,  we  shall  expand  the  total 
reflection  coefficient  In  a  series  of  partial  reflection  coefficients 
which  have  a  straightforward  physical  ray-phase  interpretation.  Having 
aone  this,  we  shall  evaluate  the  resulting  series  of  Sommerfeld  inte¬ 
grals  by  the  method  of  steepest  descents.  The  results  are  more  accurate 
than  that  of  ray  theory  and  provide  reflected  fields  readily  associated 
with  different  ray  types. 

1.3  Q/ervlew 

Chapter  2  focuses  on  the  mathematical  details  of  the  analytic 
evaluation  of  the  field  in  the  layer  using  the  L.U.A.A.  Satisfying  the 
boundary  conditions,  the  total  plane  wave  reflection  coefficient  will  be 
then  evaluated.  Finally,  the  Sommerfeld  integral  will  be  evaluated  us- 


iny  the  method  of  steepest  descents.  The  accuracy  of  these  methods  is 
also  demonstrated  by  comparison  with  numerical  results.  Two  types  of 
Interface  waves,  one  along  the  ocean  bottom  interface  and  the  other 
along  the  subbottom  Interface,  are  also  discussed.  An  equation  for  the 
interface  wave  along  the  subbottom  for  the  n-th  multiple  reflection  is 
derived. 

In  Chapter  3,  using  the  mathematical  tools  developed  in  Chapter 

2,  we  examine  a  wide  range  of  geoacoustic  models  by  parametrically  vary¬ 
ing  the  layer  thickness  and  other  geoacoustic  parameters,  as  well  as  the 
source-receiver  geometry.  We  shall  also  discuss  a  simple  but  accurate 
analytic  expression  for  the  sound  velocity  profile,  given  in  the  fbrm  of 

3 

c  -linear  in  depth,  based  on  measurements  in  order  to  make  our  geo¬ 
acoustic  models  more  realistic. 

Chapter  4  Is  devoted  to  the  discussion  of  caustics  due  to  three 
different  canonical  profile  shapes,  the  2/c  ,  c  and  c  -linear  pro¬ 
files.  Caustic  equations  are  derived  for  these  profile  shapes.  It  is 
shown  that  the  top  portion  of  the  sound  velocity  profile  In  the  layer  is 
most  critical  for  caustic  formation  based  upon  our  geoacoustic  models. 
The  fields  near  the  caustics  are  also  derived  and  examined  for  these 
profile  shapes.  It  will  be  shown  that  the  amplitude  variation  of  the 
fields  near  the  caustic  solely  depends  on  the  curvature  of  the  caustic. 

Finally,  In  Chapter  5,  real  data  collected  from  field  experiments 
will  be  analyzed  using  the  technique  developed  in  Chapter  2  and  dis¬ 
cussed  using  the  knowledge  obtained  from  the  model  studies  in  Chapter 

3.  It  will  be  demonstrated  that  the  recovery  of  geoacoustic  parameters 
from  measurements  has  a  strong  connection  to  the  geometric  orientation 


of  source  and  receiver. 
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CHAPTER  2:  Mathematical  Formulation 


l 

Analytic  solutions  to  the  problem  we  have  described  in  Chapter  1 
will  be  discussed  in  great  detail  in  this  Chapter.  Our  mathematical 
treatment  for  the  evaluation  of  the  Sommerfeld  integral  given  by  Eq. 
(1.1)  is  based  on  asymptotic  expansion  techniques  vrftich  assume  that  the 
receiver  is  many  wavelengths  away  from  the  source. 

Our  asymptotic  approximation  to  the  reflected  fields  is  twofold. 
The  first  step  is  to  compute  the  total  reflection  coefficient,  which  is 
the  weighting  factor  inside  the  Sommerfeld  integral,  by  using  Langer 
uniform  asymptotic  approximations.  The  results  are  accurate  and  yet  the 
analytic  solutions  are  general  enough  to  treat  various  typical  geoacous¬ 
tic  model  si  The  second-step  asymptotic  approximation  is  to  evaluate  the 
Sommerfeld  integral  by  the  method  of  steepest  descents  with  the  re¬ 
flection  coefficients  already  evaluated.  The  final  solutions  for  the 
reflected  fields  are  completely  analytic  and  are  functions  of  the  geo¬ 
acoustic  parameters  and  source -receiver  geometries. 

2.1  Computation  of  the  Total  Reflection  Coefficient  Using  the  Langer 
Uniform  Asymptotic  Approximation 

The  total  reflection  coefficient  for  the  proposed  geoacoustic 
model  shown  in  Fig.  1.5  will  be  given  after  general  solutions  to  the 
acoustic  velocity  potential  in  the  sediment  layer  are  evaluated  using 
the  Langer  uniform  asymptotic  approximation  (L.U.A.A.  )[10] .  Then  it 
will  be  shown  that  the  resulting  closed  form  of  the  total  reflection 


coefficient  may  be  reduced  to  ray  theory  resul  ts  for  certain  regions  of 
incident  angles  where  ray  theory  is  valid.  This  total  reflection  coef¬ 
ficient  based  on  the  L.U.A.A.  will  be  compared  to  a  numerical  result 
generated  by  the  Thoms  on- H  ask  el  1  propagator-matrix  methodCl]. 

2.1.1  General  Solutions 

We  will  refer  to  the  water  column,  the  layer  and  the  basement  as 
medium  0,  1  and  2  respectively.  First  we  consider  the  general  solutions 
in  medium  1.  The  acoustic  time  dependent  velocity  potential  W(x,y,z,t) 
satisfies  the  the  following  inhomogeneous  wave  equation 

2 

(  v2  -  - W(r,z,t)  =  0 

Cj  (z )  dt * 

After  the  harmonic  time  dependence  e_1ut  is  assumed,  W(r,z,t)  can  be 
written  as 


W(r ,z,t)  =  ( p(r,z)  e~iut 

After  the  horizontal  terms  are  separated  from  <3?(r,z ),  we  are  left  with 
the  z-dependent  potential  Fj  (z )  which  satisfies  the  onedimensional 
Helirholtz  equation. 


d2F,(z) 


A2(z)  Fj(z )  =  0 


(2.1) 


where 


u  =  angular  frequency  (rad/s) 


Q2(Z)  *-J- 


cj(z) 


(2.2) 


Here  c^(z)  is  the  sound  velocity  profile  in  the  medium  1  as  a  function 
of  depth  and  p  is  the  horizontal  slowness  or  ray  parameter.  The  hori¬ 
zontal  wave  number  kr  can  be  given  in  terms  of  p  as  kr  =  up  and  the 
vertical  wave  nuirber  k2  as  kz  =  uQ. 

Solving  for  F(z)  analytically  for  a  given  Cj(z)  is  in  general 
very  difficult  except  for  a  few  examples[5],  A  more  desirable  solution, 
which  does  not  have  to  be  exact,  but  rather  flexible  in  dealing  with 
various  acoustic  properties  in  ocean  bottom,  may  be  found  in  a  form  of 
asymptotic  approximation.  We  will  now  consider  the  following  differen¬ 
tial  equation: 


+  <Aq2(z)  -  E(z )]  G(z )  =  0 
c tr 


(2.3) 


where  E(z)  is  given  by 


E(»  .  JC  ;  f  .  ,1/6  Q-1/2 

V 

Here  n  is  known  as  the  phase  integral  and  is  given  by 


(2.4) 


n(z)  = 


£ 

j*  Q(z  ' )  dz  ' 


(2.5) 


where  z^  is  the  turning  depth  defined  by  the  condition. 
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In  general,  when  absorption  is  included  in  the  sound  velocity  pro¬ 
file,  Cj(z)  is  a  complex  function.  Thus,  zt  may  be  a  complex 
nunter[6j.  The  phase  integral  in  the  complex  z-plane  can  be  evaluated 
along  a  contour  such  as  one  shewn  in  Fig.  2.1.  The  branch  cut  due  to 
the  radical  Q  is  present  vrtiere  the  imaginary  part  of  Q  is  zero.  Since 
the  contour  always  originates  at  the  branch  point  zt,  there  is  no 
branch  cut  contribution  involved  and  the  integration  is  quite  straight¬ 
forward. 

The  basic  argument  of  the  L.U.A.A.  is  that  if  E(z)  approaches  zero 
as  <V[dc(z)/dz]  becomes  very  large  for  any  given  p,  then  G(z)  asymptoti¬ 
cally  approaches  the  desired  solution  F(z).  Indeed  it  is  possible  to 
obtain  such  an  E(z)  as  long  as  the  following  conditions  are  satisfied: 


1.  There  is  only  one  first-order  turning  point  for  a  given  p  in 
the  profile.  A  first-order  turning  point  is  defined  as  one 
for  which  Q(zt)  *  0,  but  O'  (2 1 )  /  0.  Here  the  prime 
denotes  the  derivative  with  respect  to  z. 

2.  The  profile  c^z)  must  be  smooth  [i.e.  the  condition  given 
by  Eq.  (A. 7)]  near  the  turning  point.  Since  this  condition 
must  be  met  for  any  p  and  associated  zt,  c^(z)  must  be 
smooth  everyvrfiere. 

In  brief,  c^(z)  must  be  a  smooth,  monotonically  increasing  func¬ 
tion.  The  acoustic  properties  within  a  layered  structure  of  the  ocean 
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bottom  typically  satisfy  these  conditions. 


The  exact  solutions  of  Eq.  (2.3)[9]  are  given  by 

Ga(z)  =  Q-1/2v1/4  Ai(-v)  (2.7a) 

G2(z)  =  Q-1/2v1/4  Bi(-v)  (2.7b) 

where  v  =  t|  untz)]2^3  (2.8) 

Therefore,  the  z-component  of  the  field  in  the  layer  for  the  given 
p  becomes  asymptotically 

Fj ( z  )  “  Q"1/2v1/4[B1Ai(-v)  ♦  B2  Bi(-v)]  as  u  >  co  (2.9) 

The  error  analysis  of  the  L.U.A.A.  is  discussed  in  Appendix  A. 

Since  medium  0  is  assumed  to  be  homogeneous,  the  z-component  of 
the  general  solutions  for  the  velocity  potential  can  be  easily  obtained 
ana  is  given  by 

Fq(z)  =  Ax  eiuYoz  +  Ag  e'7uY0z  (2.10) 

where  y  =  (-L-  -  p2)1/2  (2.11) 

c0 

Similarly,  in  medium  2 

F2(z)  *  Cj  e'juY2z  (2.12) 
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where 


(2.13) 


F 


i 


I 

I 

l 


y 


.  p2i1/2 


Note  that  there  is  only  one  term  in  Eq.  (2.12)  in  order  to  satisfy  the 
radiation  condition.  Aj,  Ag,  Bj,  B2  and  Cj  are  arbitrary  constant  co¬ 
efficients  for  a  given  p  and  remain  to  be  found  from  appropriate  bounda¬ 
ry  conditions,  which  will  be  discussed  shortly.  After  multiplication  by 
the  horizontal  component,  the  potential  is  given  by 


<£t(r,  z)  =  F.(z)  e1tjpr;  i  =  0,  1,  2  (2.14) 

in  the  i-th  medium  for  a  given  p. 

2.1.2  Boundary  Conditions 

Since  the  pressure  and  vertical  component  of  the  motion  of  the 
particles  must  be  continuous  across  the  boundaries,  p^  and  d^/dz  must 
be  equal  at  the  interfaces  in  order  to  satisfy  these  conditions  respec¬ 
tively.  Thus  we  can  obtain  the  following  equations: 

At  z  =0  (the  ocean  bottom) 

pq(A1  ♦  A2)  ==  pjCBjGjtO)  ♦  B2G2(0)]  (2.15) 

iuVAl  ‘  A2}  s  BiGi (0)  ♦  B2G^(0)  (2.16) 

At  z  =  H  (the  subbottom) 

p1CB1G1(H)  +  B2  Gg  ( H )  3  =  P2C1e-iur2H  (2.17) 
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.  ft 


p 


BjG^H)  ♦  B2GJ(H)  =  -  1CiUy2  e'1'UY2H 


(2.18) 


where 


G  x(z  >  =  Q‘1/2v1/4Ai(-v) 


G2(z)  =  Q‘1/2v1/4Bi(-v) 


(2.19a) 

(2.19b) 


Solving  for  Aj  assuming  a  unit  amplitude  in  the  incident  wave  (A,,  =  1), 
we  finally  obtain  the  total  reflection  coefficient  for  this  geoacoustic 
ocean  bottom  model , 


R(p)  = 


(A, 


(2) 


<*«» 


(A»>  -B «»,,*«« 


(*J2)  -s'2))  (a'J1  -B‘»)  -  .Bj1’)  (a‘2) 


-b'2)) 


(2.20) 


where 

Aq2)*  (2^  *  Tp1y0W(1)*{2}  (z  ac  0) 

( 2.21a  ,b) 

B<1,*(2)  a-P0Q0v-1/2W(1)*(2)  (z  *  0) 

( 2.22a  ,b) 

A<1J‘(2)  *  iPlY2W(1)’(2)  (z  *  H) 

(2.23a  ,b) 

BH1),<2)  -p^1/2w,(1),(2)  (*■»> 

(2.24a  ,b) 

wi  th 

Q0  =  Q(z=0),  Qh  =  O(z-H) 

Vq  m  v(ZsO),  V^|  =  v(ZsH) 

w(1)’(2)(v(z))  *  Ai(-v)  ~  iBi(-v) 

(2. 25a, b ) 

H.(  D,( 2)  (- (T ))  (dAi(-v)  -  dBi(-v)j 

Oy  dv 

(2.26a,b ) 
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Here  are  known  as  Fock  funct1ons[9].  The  solutions  may  be 

given  In  terms  of  Hankel  functions  of  the  l/3rd  order.  However,  the 
advantage  of  using  Fock  functions  Is  that  because  they  are  linear  com¬ 
binations  of  Airy  functions,  it  Is  easier  to  evaluate  them  numerically 
than  Hankel  functions.  In  addition,  similar  to  Hankel  functions,  their 
asymptotic  expansions  (see  Appendix  C)  behave  like  propagating  plane 
waves.  In  general,  and  are  not  complex  conjugates  of  one 
another  since  v  as  well  as  Ai  and  Bi  may  be  complex  nunfoers. 

This  asymptotic  approximation  of  the  total  plane  wave  reflection 
coefficient  given  by  Eq.  (2.20)  is  valid  for  any  shape  of  the  sound 
velocity  profile  in  the  sediment  layer  as  long  as  it  is  a  slowly  and 
monotonically  Increasing  profile.  Not  only  these  conditions  are  met  in 
geoacoustic  modeling  of  a  typical  ocean  bottom,  but  also  this  new  plane 
wave  reflection  coefficient  allows  us  to  do  more  detailed  geoacoustic 
modeling  by  providing  some  freedom  in  the  choice  of  profile  shapes.  The 
total  plane  wave  reflection  coefficient  based  on  this  formulation  Is 
compared  to  the  numerical  result[13]  generated  by  the  Thomson-Haskel  1 
Propagator-tatrix  method[l]  for  a  typical  case  of  Interest  (cf.  Fig. 

2.2)  for  a  frequency  of  220  Hz.  As  shown  in  Fig.  2.3,  their  agreement 
is  so  good  that  one  may  not  distinguish  one  from  the  other. 
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Fig.  2.3  The  total  reflected  coefficients  for  the  model  shown  in  Fig. 

2.2.  (a)  The  Thomson- H  ask  el  propagator-matrix  method  with  a 

step  size  of  Z-m  in  the  layer,  (b)  The  Langer  uniform 
asymptotic  approximation. 


2.2  Physical  Properties  of  Reflection  Coefficients 


2.2.1  Expansion  of  the  Total  Reflection  Coefficient 
In  this  section,  the  total  reflection  coefficient  obtained  using 
the  L.U.A.A.  will  be  decomposed  into  an  infinite  sum  of  partial  reflec¬ 
tion  coefficients.  It  will  be  shown  that  each  term  represents  its  own 
unique  physical  property  analogous  to  ray  theory  for  special  cases  where 
the  turning  depths  are  away  from  both  of  the  interfaces. 

The  procedure  of  the  expansion  is  analogous  to  a  ray  expansiontll ] 
which  decomposes  the  total  reflection  coefficient  into  an  infinite  sum 
of  reflection  coefficients:  First,  we  divide  the  numerator  and  denomi¬ 
nator  of  Eq.  (2.20)  by 

(A<*>  *8<2>Ka'1> 


to  obtain 


R(p )  = 


a(2)  d  (2) 
A0  "  B0 
a(2)  7  „(2) 
Ao  Bo 


(A(1)  _  b(1)){a(2) 
1*0  -  Bo 


(A*21  *b‘2,)(a‘11 


77^ 


,A(1)  .  b(1)|(A(2> 
°o  h*h 


-  bh2)  > 


-1 


<A<«  +  B0Z,)(AHlJ  •  BH1}) 


(2.27) 


Second,  we  perform  a  binomial  expansion  of  the  second  bracket. 
This  process  may  be  justified  by  introducing  a  small  positive  imaginary 
part  in  c^(z)  which  makes  the  absolute  value  of  the  second  term  always 


smaller  than  unity,  and  letting  the  imaginary  part  go  to  zero  after  the 
expan-  sion.  After  a  few  more  steps,  R(p)  becomes 

00 

R(p)  =  Rq  Rn  (2,28) 
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In  order  to  obtain  a  physical  interpretation  of  each  term,  we  sub¬ 
stitute  the  asymptotic  forms  of  the  Airy  functions.  The  conditions  for 
which  the  Airy  functions  may  be  expanded  are  equivalent  to  the  condition 
that  the  turning  depth  for  a  given  horizontal  slowness  p  is  many  wave¬ 
lengths  away  from  either  boundary.  The  first  term,  Rq  can  then  be  re¬ 
duced  to 
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(2.31) 


This  is  simply  the  Rayleigh  plane  wave  reflection  coefficient  in  medium  0 
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over  the  homogeneous  half  space  vrfiich  is  characterized  by  c^(z=0)  and 
This  is  the  plane  wave  reflection  coefficient  one  would  expect  to 
obtain  for  the  specularly  reflected  ray  using  ray  theory.  Similarly, 
for  large  v,  R^  can  be  reduced  to 


T01I10R12e21“','0"'H’ 


t  t  P2l<Jnn  "  1ir//2 
‘oi 1 10e  u 


2t<  H 


°  >  zt  >  H 


(2.32) 


where  TQ^,  T1q  and  R12  are  local  plane  wave  transmission  coefficients 
and  reflection  coefficient  across  the  medium  0  to  1 ,  1  to  0  and  1  to  2, 
respectively.  This  is  the  result  one  would  obtain  from  ray  theory  for  a 
ray  which  reflects  off  the  subbottom  or  turns  within  the  layer  once,  de¬ 
pending  on  the  location  of  the  turning  depth.  Notice  the  -ir/2-phase 
shift  associated  with  the  refracted  ray,  a  term  which  is  often  introduc¬ 
ed  as  a  correction  in  pure  ray  theory.  In  exactly  the  same  manner  Rn 
can  be  reduced  to  the  following  asymptotic  form  for  n  >  1 . 
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vrfiere  R's  and  T's  with  two  integer  subscripts  are  local  plane  wave  re¬ 
flection  and  transmission  coefficients  at  the  interface  indicated  by 
their  subscripts.  From  the  phase  in  Eq.  (2.33),  it  suggests  that  Rp 
behaves  asymptotically  in  n  like  the  n-th  multiple  which  reflects  or  re¬ 
fracts  n  times  in  the  layer.  Rj  is  a  special  case  of  Rn  with  n  =  1. 


because  of  the  physical  interpretation  obtained  by  using  the  asymptotic 
approximation,  we  may  call  Rq,  and  Rn  the  specular,  the  primary  and 
the  n-th  multiple,  respectively.  The  schematics  of  their  physical  prop¬ 
erties  are  shown  in  Fig.  2.4. 

We  note  that  while  reflection  coefficients  based  on  ray  theory 
fail  when  the  turning  depth  is  very  close  to  either  interface,  these  new 
reflection  coefficients  based  on  the  L.U.A.A.  are  uniformly  valid  for 
all  values  of  p  regardless  of  the  locations  of  the  corresponding  turning 
aepths.  This  is  one  of  the  most  powerful  features  of  the  Langer  reflec¬ 
tion  coefficients. 


2.2.2  Comparison  with  Ray  Theory  Results 
Interesting  features  of  the  Langer  reflection  coefficients  as 
compared  to  reflection  coefficients  based  on  ray  theory  will  be  dis¬ 
cussed  next.  In  Figs.  2. 5-2. 7,  the  specular,  primary  and  second  multi¬ 
ple  Langer  reflection  coefficients  are  compared  to  the  ray  theory  re¬ 
sults  for  a  geoacoustic  model  shown  in  Fig.  2.2  at  a  frequency  of  220  Hz. 

Because  of  the  small  sound  velocity  drop  at  the  ocean  bottom,  we 
observe  a  minimum  amplitude  in  the  specular  reflection(cf.  Fig.  2.5a)  at 
a  particular  incident  angle.  According  to  ray  theory  we  would  expect 
the  specular  reflection  coefficient  to  be  zero  at  this  angle(cf.  Fig. 
2.5b),  which  is  the  Brewster  angle  (angle  of  intromission)  associated 
with  two  homogeneous  media  where  the  upper  half-space  has  the  properties 
p0  and  cQ,  and  the  lower  half-space  has  the  properties  >  pQ  and 
Cj(z  *  0)  <  Cq  [4].  However,  the  specular  reflection  coefficient  does 
not  become  identically  zero  according  to  the  L.U.A.A.,  but  only  becomes  a 
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Fig.  2.4  Typical  reflecting  rays  in  the  bottom. 
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5  The  specular  reflection  coefficient.  A  comparison  between 
(a)the  L.  U.  A.  A.  and  (b)ray  theory  for  the  model  shown  in 
Fig.  2.2.  The  frequency  is  220  Hz. 
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Fig.  2.7  The  second  multiple  reflection  coefficient.  A  comparison 
between  (a)the  l.  U.  A.  A.  and  (b)ray  theory  for  the  model 
shown  in  Fig.  2.2.  The  frequency  Is  220  Hz. 


minimum.  This  Is  due  to  the  fact  that  the  positive  sound  velocity 
graaient  in  the  sediment  layer  prohibits  total  wave  energy  transmission 
even  at  the  Brewster  angle.  This  discrepancy  is  also  observed  in  the 
plane  wave  reflection  coefficients  for  the  mul tipi es (cf.  Fig.  2.7). 

Another  Interesting  feature  of  the  Langer  reflection  coefficients 
can  be  seen  for  the  geoacoustic  model  shown  in  Fig.  2.8.  Here  the  sound 
velocity  is  continuous  at  the  water-bottom  interface,  but  there  is  still 
a  velocity  graaient  in  the  bottom  and  the  densities  in  both  media  are 
the  same.  Ray  theory  predicts  no  specular  reflection  or  multiple  re¬ 
flections  associated  with  the  water-bottom  interface  because  there  is  no 
acoustic  impedance  discontlnui ty  exists  across  that  Interface. 

However,  according  to  the  Langer  reflection  coefficients,  as  shown 
in  Fig.  2.9,  one  may  still  obtain  non-zero  reflection  coefficients  near 
grazing  incidence.  These  unexpected  results  may  be  explained  by  the 
fact  that  the  Langer  reflection  coefficients  take  Into  account  not  only 
the  discontinuity  in  the  acoustic  impedance  but  also  the  discontinuity 
In  the  derivative  of  the  acoustic  impedance  with  respect  to  z.  This  is 
another  example  where  the  Langer  reflection  coefficients  demonstrate 
their  improved  accuracy  over  ray  theory  based  reflection  coefficients 
for  a  stratified  slowly  varying  medium. 

2.3  Asymptotic  Expansion  of  the  Sommerfeld  Integral  by  the  Method  of 

Steepest  Descents 

Our  goal  Is  to  compute  the  reflected  fields  due  to  a  harmonic 
point  source,  namely  to  obtain  the  leading  order  asymptotic  behavior  of 
the  Sommerfeld  integral,  Eq.  (1.1),  by  means  of  as>mptotic  approxima¬ 
tion.  So  far  we  have  found  the  closed  form  of  R(p)  in  Eq.  (1.1)  for  our 


geoacoustic  models,  and  further  how  to  break  up  R(p)  into  an  infinite 
number  of  terms  each  of  which  has  a  physical  interpretation  in  terms  of 
ray  theory. 

After  substituting  Eq.  (2.28)  in  Eq.  (1.1),  we  obtain  a  series  of 
Sommerfeld  integrals  each  of  which  contains  Rn<  An  asymptotic  form  of 
each  Sommerfela  integral  can  provide  us  an  analytic  expression  for  the 
reflected  field  as  a  function  of  horizontal  range  coupled  with  the  phys¬ 
ical  properties  of  Rj,  associated  wi th  the  integrand,  khile  the  co¬ 
herent  sum  of  these  asymptotic  forms  of  the  Integrals  provides  the  total 
reflected  field,  we  can  also  evaluate  how  the  individual  reflected 
fields  contribute  to  the  total  reflected  field  along  horizontal  range 
for  various  sets  of  geoacoustic  parameters  and  source-receiver  geometri¬ 
es  . 

Thus,  after  the  Hank  el  function  is  replaced  with  its  asymptotic 

form. 
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from  which  we  obtain  the  reflected  field  representing  the  physical  prop- 
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erty  that  Rr  carries. 


2.3.1  Mathematical  Treatment  before  the  Method  of  Steepest 
Descents  Is  Applied 

Fig.  (2.10)  shows  the  complex  p-plane  with  branch  points  at  p  = 
l/cQ  and  l/c2<  (Only  the  l/cQ-branch  point  is  present  for  PQ  contain¬ 
ing  Rq. )  The  branch  cuts,  known  as  the  E.J.P.  cuts[8] ,  satisfy  the 
conoition  that  the  imaginary  parts  of  the  radicals,  p,  and  y2  are 
zero.  Notice  that  they  are  the  only  radicals  that  we  must  be  concern¬ 
ed  with  in  the  process  of  deforming  the  integration  contour  to  eval¬ 
uate  Eq.  (2.34)  even  after  the  expansion  of  the  total  reflection  coef¬ 
ficient.  Unlike  a  slab  configuration  of  homogeneous  media[10] ,  i.e.  a 
homogeneous  layer  between  two  different  homogeneous  half  spaces,  we  do 
not  create  additional  branch  points  as  a  result  of  the  ray  expansion  of 
the  total  reflection  coefficient.  Qq  and  Qh  may  appear  as  radicals  in 
Rn  at  first  glance  but  with  v-1^2  there  is  no  antoiguity  in  sign. 

We  shall  apply  the  method  of  the  steepest  descents[2]  to  eval¬ 
uate  P„.  In  order  to  do  so  we  must  rewrite  P„  in  the  form  of 
n  n 


I„  -  f  Fn(p)  eMn(P}  dp  ;  u  >  co  (2.35) 

•'-CD 

where  F(p)  must  be  a  slowly  varying  function  of  p,  compared  to  the  com¬ 
plex  exponential  term  in  order  for  the  method  to  be  valid.  When  n  4  0, 
this  would  not  be  the  case  if  we  made  the  direct  substitutions 


Fn(p)  =  R 


0n(p)  =  Y0( Z+ZQ)  ♦  pr 
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because  Rn  is  in  general  a  very  fast  oscillatory  function  (See  Fig. 

2.6  or  Fig.  2.7  for  example).  In  order  to  convert  Eq.  (2.34)  into  the 
form  of  Eq.  (2.35)  satisfying  the  above  requirement,  the  following  sub¬ 
stitutions  are  suggested: 

(2.36) 

VP)  =  Vz  +  zo)  +  Pr  +  2nw  (2.37) 

where 
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0  >  zt  >  H  (2.39a) 

z  t  <  H  (2.39b) 


in  trfiich  the  phase  term  of  Rn  is  estimated  from  the  asymptotic  behav¬ 
ior  of  Rn,  which  is  identical  to  the  ray  phase,  then  separated  from 
R,,  and  added  to  the  complex  exponential  term  in  the  integrand.  Fig. 

2.11  shews  Vn  vs.  incident  angle  for  the  primary  reflection  coeffi¬ 
cient  of  the  geoacoustic  model  shown  in  Fig.  2.2.  Compared  to  Fig.  2.6, 
the  phase  curve  has  been  significantly  stabilized.  The  only  appreciable 
phase  variation  occurs  between  critical  incidence  and  the  angle  at  which 
the  ray  completely  turns  in  the  layer.  Since  Vn  approaches  Rp  as  u 
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11  The  primary  reflection  coefficient  for  the  model  shown  in 
Fig.  2.2  after  the  phase  term2nn  is  removed. 
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becomes  very  smalltcf.  Eq.  (2.38)],  the  substitutions  given  by  Eq. 
(2.36)  through  Eq.  (2.38)  are  valid  even  when  nun  is  not  large. 


2.3.2  Evaluation  of  Pn  by  the  Modified  Method  of 

Steepest  Descents  -  Saddle  Point  Contribution 

Strictly  speaking,  a  direct  application  of  the  method  of  steepest 
descents  to  the  Sommerfeld  integral  with  the  Langer  reflection  coeffi¬ 
cient,  which  is  a  function  of  p  and  u ,  is  not  valid.  However,  by  ex¬ 
cluding  the  second  derivative  of  Vn  in  the  expansion,  we  may  still 
apply  the  basic  technique  to  obtain  analytic  solutions  to  the  Sommerfeld 
integral  whose  accuracy  is  improved  over  the  results  obtained  by  ray 
theory.  The  detailed  derivations  and  justifications  are  discussed  in 
Appendix  B,  ana  here  we  only  introduce  the  necessary  formulation  to 
obtain  the  asymptotic  solution  to  the  Sommerfeld  integral  of  the  n-th 
term. 

First,  we  must  find  a  saddle  point  ps  which  satisfies 
Ps 

(ps)  =  -  ~  (z  +  zQ)  ♦  r  ♦  2nu'  (p$)  =  0  (2.40) 

This  is  equivalent  to  locating  an  eigenray  connecting  two  fixed  points. 

In  some  cases,  there  could  be  well  separated  multiple  saddle  points  (but 
not  a  caustic),  found  for  the  same  source-receiver  geometry. 

Second,  we  deform  the  integration  contour  passing  through  the  sad- 
ale  point(s)  along  which  the  real  part  of  *5  is  constant  and  the  imagi¬ 
nary  part  is  positive.  Such  a  deformed  contour  is  shown  in  Fig.  2.12. 

As  a  result,  the  complex  exponential  term  now  becomes  a  decaying  term 
along  this  deformed  contour  with  a  maximum  at  the  saddle  point.  There- 
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fore.  if  the  amplitude  of  the  integrand  is  a  well  behaved  function,  the 
integration  may  be  approximated  by  integrating  through  a  small  domain 
near  the  saddle  point  (marked  with  ")("  in  Fig.  2.12),  and  the  rest  of 
the  contribution  asymptotically  becomes  zero. 

Next  we  expand  the  integrand  in  a  Taylor  expansion  in  p  about  ps> 
and  Eq.  (2.35)  becomes 


P„  =  f  dp  [Fn(ps)  •  F'n(ps)(p  -  ps)  -  i  F„- {ps > (p  -  p$)2] 

-  CO 
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2  3  •  4. 

ana  the  asymptotic  expansion  of  Eq.  (2.34)  due  to  the  saddle  point  Pns 
up  to  the  second  order  term  becomes[2] 
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The  above  formula,  Eqs.  (2.41)  through  (2.55),  are  also  valid  for 
complex  sound  veloci  ties  to  include  absorption.  In  this  case,  a  complex 
saddle  point  must  be  found  from  Eq.  (2.40)  where  the  phase  integral  u  is 
evaluated  for  a  complex  turning  depth  zt.  We  also  point  out  that  Eq. 
(2.41)  is  the  contribution  from  a  single  saddle  point.  Thus,  i  f  more 
than  a  single  saddle  point  exists,  their  contributions  must  be  evaluated 
and  summed. 

If  N  is  set  equal  to  zero[cf.  Eq.  (2.42)],  the  r es ul t  b ecomes  a 
ray  solution  except  Rn  is  more  accurately  evaluated  by  the  L.U.A.A. 

It  is  necessary  to  include  the  correction  term  N  in  order  to  be  con¬ 
sistent  with  the  order  of  accuracy  provided  by  the  Langer  formulation 
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which  is  valid  up  to  order  1/u  (See  Appendix  A).  However,  using  the 
same  argument,  evaluation  of  higher  order  terms  in  Eq.  (2.42)  or  more 
accurate  integration  of  Eq.  (2.34)  by  a  numerical  method  would  not  in¬ 
crease  the  accuracy  of  the  result  given  by  Eq.  (2.42). 

2.3.3  Interface  Waves 

Dealing  with  two  different  types  of  interfaces  in  our  geoacoustic 
models,  we  consider  two  physically  and  mathematically  different  types  of 
Interface  waves;  one  propagates  along  the  ocean  bottom  interface  where 
the  medium  containing  the  incident  wave  is  homogeneous  and  overlies  an 
inhomogeneous  halfspace,  while  the  other  propagates  along  the  ocean  sub¬ 
bottom  where  the  medium  containing  the  incident  wave  is  inhomogeneous 
ana  overlies  a  homogeneous  halfspace. 

The  first  type  of  interface  waves,  known  as  an  interference  head 
wave,  has  been  extensively  studied  by  Cerveny  and  Ravindra[5].  It  is 
the  result  of  a  coherent  sum  of  many  multiples  which  have  closely  lo¬ 
cated  shallow  turning  depths  (cf.  Fig.  2.13). 

By  defini tion[7] ,  the  interference  head  wave  propagates  along  the 
interface  such  that  the  upper  halfspace  is  homogeneous  with  the  sound 
velocity  cQ  and  the  lower  halfspace  has  the  acoustic  property  as 
c j^z)  =  a  +  bz  where  a  >  cQ  and  b,  z  >0.  However,  unless  a  limit¬ 
ing  case  where  b  approaches  0  in  Cj(z)  is  considered,  the  mathematical 
arguments  are  Identical  for  the  case  where  c^(0)  Is  greater  or  less 
than  cQ,  and  both  are  the  contributions  from  the  multiples. 

We  have  observed  that  there  are  no  additional  branch  points  cre¬ 
ated  as  the  result  of  the  expansion  of  the  total  reflection  coeffi- 
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cient.  Thus,  unlike  head  waves  from  the  branch  cut  contribution^]  for 
two  homogeneous  halfspaces,  the  Interference  head  waves  are  special 
cases  of  the  saddle  point  contributions  to  the  total  reflected  field. 
Therefore,  the  amplitudes  of  the  interference  head  waves  are  proportion¬ 
al  to  the  reciprocal  of  the  range,  not  the  range  squared. 

The  presence  of  the  second  type  of  interface  wave  Is  due  to  the 
branch  cut  contribution  as  the  S.D.P.  crosses  the  branch  point  at  p  = 
l/c2  (but  never  at  p  =  1/ Cq),  and  the  contour  must  be  deformed  to  go 
arouna  this  branch  cut  so  that  the  end  points  of  the  integration  path 
may  remain  on  the  correct  Riemann  sheettcf.  Fig.  2.14).  The  technique 
to  evaluate  the  asymptotic  behavior  of  the  contour  integration  around 
the  branch  cut  is  completely  analogous  to  that  of  Lanb  's  problem[l], 

[12 j  where  the  method  of  steepest  descents  is  applied.  The  asymptotic 
approximation  of  the  interface  wave  due  to  the  n-th  multiple  including 
the  primary  reflection,  PnL  is  given,  as  expected,  by 
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where  L  Is  the  horizontal  distance  the  interface  wave  propagates  along 
the  subbottom  for  each  subbottom  interaction (cf.  Fig.  2.15).  As  r  be¬ 
comes  large  so  does  nl,  and  PnL  algebraically  decays  proportional  to 
l/r*.  Although  Eq.  (2.56)  is  the  leading  order  term  of  the  asymptotic 
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Fig.  2.15  Ray  interpretation  of  the  interface  waves  propagating  along 
the  subbottom  interface,  (a)  The  primary  reflection,  (b )  the 
2nd  multiple  reflection  and  (c)  the  3rd  multiple  reflection. 


expansion  of  PnL  •  U  is  of  order  1/u.  Thus,  this  is  all  we  can  expect 

from  the  Langer  reflection  coefficients. 

The  total  subbottom  interface  wave  contribution  is  given  by  the 

sum  of  the  branch  cut  contributions  of  P„.  as  their  S.D.P.'s  cross  the 

nL 

branch  point.  Again,  like  the  saddle  point  contribution,  Eq.  (2.56)  is 
valid  for  any  inhomogeneous  layer  as  long  as  it  satisfies  the  conditions 
necessary  for  the  Langer  reflection  coefficient. 

2.4  Limitation  of  the  Asymptotic  Expansion  of  Pn 

The  total  reflected  fields  at  two  different  values  of  z  +  Zq  are 
computed  from  the  geoacoustic  model  shown  in  Fig.  2.7.  The  results  are 
compared  to  the  ones  generated  by  a  numerical  Hankel  transform  method 
L12](cf.  Fig.  2.16).  They  show  good  agreement  except  at  a  few  rigions 
of  horizontal  ranges  where  the  asymptotic  solutions  to  the  primary  and 
multiple  reflections  break  down  at  the  critical  incidence  and  caustics. 

As  discussed  in  Secs.  2.3.1  and  2.3.2,  in  order  for  the  method  of 
steepest  descents  to  be  valid,  F(p)  must  be  a  smooth  well-behaved  func¬ 
tion.  Although  F(p)  can  be  made  sufficiently  smooth  at  most  incident 
angles  by  separating  out  the  complex  phase  term,  the  rapid  change  in 
amplitude  and  phase  near  critical  incidence  is  still  a  problem.  Since 
V'n(p)  at  critical  incidence  becomes  large,  the  correction  term  N  be¬ 
comes  unstable.  This  is  also  true  for  reflections  at  very  small  graz¬ 
ing  angles.  The  second  problem  is  that  N  blows  up  as  0''(ps)  ap¬ 
proaches  zero.  This  is  the  condition  corresponding  to  caustic  formation. 

In  the  complex  p-plane,  the  first  problem  is  due  to  the  fact  that 
the  saddle  point  is  located  too  close  to  the  branch  point,  1/^  or  1/Cq. 
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A  comparison  between  the  asymptotic  method  and  the 
numerical  Hankel  transform  method  at  z  +  Zg  =  250  i 
for  the  geoacoustic  model  shown  in  Fig.  2.Z. 


RELATIVE  AMPLITUDE  (dB) 


Fig. 


2.16b  A  comparison  between  the  asymptotic  method  and  the 
numerical  Hankel  transform  method  at  z  ♦  zg  =  150  n 
for  the  geoacoustic  model  shown  in  Fig.  2.2. 
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For  the  second  problem,  the  S.O.P.  passing  through  the  saddle  point  is 
no  longer  a  straight  line  near  the  saddle  point,  but  branches  out  with 
w/3  angles  (cf.  Fig.  2.17).  Thus  the  expansion  along  the  straight 
11ne(cf.  Fig.  2.12),  which  resulted  in  Eq.  (2.41),  is  no  longer  valid  at 
the  caustic  saddle  point  and  breaks  down. 

For  either  case  another  alternative  must  be  taken  in  order  to 
evaluate  Pn  for  such  regions  of  p.  A  uniform  asymptotic  expansion  of 
the  integral  across  the  critical  points  [3]  may  be  applied  for  these 
cases,  but  it  is  very  difficult  to  obtain  any  physical  interpretation 
from  the  results  which  involve  special  functions.  Thus,  the  argument 
for  obtaining  analytic  solutions  for  these  special  conditions  becomes 
very  weak.  Therefore,  our  choice  here  is  to  numerically  integrate 
along  a  slightly  modified  S.D.P.  as  shown  in  Fig.  2.18  under  these  cir¬ 


cumstances  . 
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Chapter  3:  Geoacoustic  Models  and  Reflected  Fields 


In  Sec.  2.1  the  total  plane  wave  reflection  coefficients  for  our 
geoacoustic  models  were  derived  using  the  Langer  Uniform  Asymptotic 
Approximation  (L.U.A.A. ).  Then,  using  the  Langer  reflection  coeffi¬ 
cients,  we  derived  reflected  fields  due  to  a  time  harmonic  point  source 
in  Sec.  2.3.  We  shall  next  derive  the  phase  integrals  for  the  profile 
shapes  of  interest. 

After  examining  the  full  analytic  solutions  for  the  reflected 
fields,  particularly  the  asymptotic  behavior  of  the  fields  at  long 
ranges,  we  can  estimate  important  parameters  to  study  various  cases.  We 
shall  then  apply  the  above  results  to  study  energy  partitioning  of 
eigenrays  for  the  ocean  bottom  geoacoustic  models. 

3.1.  Profile  Types  in  the  Sediment  Layer 

The  formulae  for  the  reflection  coefficient  given  by  Eq.  (2.20) 
and  the  asymptotic  expansion  of  the  reflected  fields  given  by  Eq.  (2.42) 
through  (2.59)  are  all  based  on  an  ocean  bottom  geoacoustic  model  con¬ 
sisting  of  a  layer  with  a  slowly  varying  sound  velocity  (of  any  shape) 
above  a  homogeneous  infinite  halfspace.  Therefore,  the  only  parameter 
vWiich  must  be  found  is  the  phase  integral  n  defined  by  Eq.  (2.5).  This 
quantity  may  be  calculated  numerically,  but  it  is  instructive  to  derive, 
where  possible,  an  analytic  expression  in  order  to  obtain  full  analytic 
solutions  for  the  reflected  fields  and  to  investigate  the  physical  ef¬ 
fects  of  the  geoacoustic  parameters  upon  the  reflected  fields.  We  shall 


limit  our  cases  to  three  different  profile  shapes  for  which  analytic 
expressions  for  the  phase  Integral  and  its  derivatives  are  obtainable. 


a)  1/Cj -Linear  Profile 
For  this  case 


Cj(z) 


-  ai  2  ’  bi 


(3.1) 


where  a^  and  b^  are  arbitrary  parameters.  This  profile  type  is  one 
of  the  most  commonly  used  for  modeling  slowly  varying  inhomogeneous 
media  because  exact  analytic  solutions  for  the  field  are  available  in 
the  form  of  Airy  functions.  Also,  the  Langer  reflection  coefficients 
become  exact  for  this  profile  type  as  £(z)  in  (2.3)  vanishes  every¬ 
where.  Even  though  one  might  argue  that  the  shape  of  this  profile  does 
not  necessarily  represent  actual  profile  shapes  observed  in  real  ocean 
bottom  sediments,  it  will  be  shown  in  Sec.  3.2.3  that  this  is  a  good 
approximation  to  the  correct  ones  as  long  as  the  layer  thickness  and 
souna  velocity  gradient  are  small. 

The  analytic  expressions,  for  the  first,  second,  third  and  fourth- 
-order  derivatives  of  eta  with  respect  to  p  for  this  profile  type  are 
given  below: 


(3.2) 
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where  Q  a  [l/c2(z)  -  p2]1^2. 


b)  Cj-1  inear  Profil  e 
For  this  case 

cx(z)  *  a2  z  +  b2 


(3.4) 


(3.5) 


(3.6) 


(3.7) 


where  ^  and  b2  are  arbi  trary  parameters. 

This  is  another  common  profile  shape,  a  straight  line  between  two 
fixed  endpoints.  Formulae  for  the  phase  integral  and  derivatives  for 
tins  profile  are  given  by 


,  1  ♦  c,(z)Q 

"2  “  aj  tcl (z ><1  -  ln(  pCjU)  i] 

dtTg  c  x  (z)Q 
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(3.11) 


(3.12) 


c)  c^-L.1near  Profile 


For  this  case 


CjU)  •  aj2.  b3 


(3.13) 


where  a-j  ana  b^  are  arbitrary  parameters. 

This  rather  unusual  linear  expression  for  the  sound  velocity  is 
added  to  our  analysis  in  order  to  investigate  more  realistic  geoacoustic 
models  derived  from  measurements^] .  The  sound  velocity  in  the  bottom 
may  be  given  as  a  function  of  one-way  travel  time  t  from  direct 
measurements  and  is  often  expressed  in  terms  of  polynomials  in  t.  From 
these  expres-  sions,  Hamilton  has  studied  sonobuoy  bottom  reflection 
data  ana  statis-  tically  obtained  a  form  for  the  sound  velocity  as  a 
function  of  depth  which  he  fitted  with  third-order  polynomial s[2].  As 
seen  in  Fig.  D.l ,  the  major  characteristic  of  the  "Hamilton  profile"  is 
that  the  sound  velocity  gradient  near  the  ocean  bottom  is  larger  than 
that  of  the  deep-  er  portion  of  the  profile  -  this  is  just  the  opposite 
of  the  behavior  for  the  case  of  the  1/c^  (z)-l  Inear  profile. 

It  has  been  found  by  the  author  that  the  "Hamilton  profile"  may  be 
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l 


well  approximated  by  a  simple  linear  expression  as  Introduced  here. 
From  the  given  coefficients  of  a  third-order  polynomial  in  t,  one  may 
obtain  a^  and  bj  in  (3.13)  from  the  formula  provided  in  Appendix  D. 
Their  agreement  Is  quite  good  even  In  the  deep  portion  of  the  profile 
where  the  maximum  error  is  expected  (See  Fig.  D.l).  The  validity  of 
this  rather  simple  approximation  is  also  discussed  in  Appendix  D. 

The  analytic  expressions  for  n  and  its  derivatives  may  be  easily 
obtained  and  are  given  by 


n3  =  - 
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Similarly,  a  profile  whose  form  is  c  -linear  wl  th  M  =  2,  4  or  5 
also  provides  simple  analytic  expression  for  the  phase  integral.  The 
general  rule  is  that  the  greater  M  is,  the  greater  the  gradient  near  the 

3 

ocean  bottom.  These  profiles  may  be  used  instead  of  the  c  -linear 


form,  depending  on  the  size  of  the  curvature  of  the  profile  shape  to  be 
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modeled.  We  choose  to  study  the  c  -linear  profile  because  it  demon¬ 
strates  the  best  agreement  with  the  measurements  reported  by  Hamilton. 

In  order  to  study  the  effects  due  solely  to  the  profile  shapes,  we 
fix  the  end  points  of  the  profiles  so  that  the  difference  in  character¬ 
istic  acoustic  impedance  across  each  interface  remains  the  same.  Once 
the  endpoints  of  the  sound  velocity  in  the  layer  are  fixed,  the  coef¬ 
ficients  of  the  linear  equation  for  each  profile  shape  are  automatically 
fixed. 

3.2  Important  geoacoustic  parameters  in  modeling 

In  this  section,  applying  the  mathematical  results  obtained  in 
Chapter  2  and  profile  shapes  chosen  for  the  layer  in  Sec.  3.1,  we  shall 
discuss  the  effects  of  geoacoustic  parameters  on  the  reflected  fields 
ano  thus  establish  guidelines  to  choose  appropriate  examples  of  geo¬ 
acoustic  models.  Since  we  have  obtained  closed  forms  for  the  reflected 
fields  based  on  the  geoacoustic  models,  we  can  discuss  the  direct  para¬ 
metric  dependence  of  our  solutions  for  some  limiting  cases. 

The  major  quantities  which  are  primarily  responsible  for  the  field 
amplitudes  for  a  given  geoacoustic  model  are  b"_1,/^(ps)  and  Rn(Ps)  Ccf. 
Eq.  (2.42)].  The  first  one  controls  the  geometric  spreading  and  the 
other  the  amplitude  of  the  plane  wave  corresponding  to  the  particular 
saddle  point.  Since  we  have  fixed  the  endpoints  of  the  layer  profile  in 
modeling,  the  boundary  conditions  at  both  of  the  interfaces  remain  the 
same.  This  leads  to  the  condition  that  the  asymptotic  behavior  of  the 
amplituoe  of  the  reflection  coefficient  Rn  is  independent  of  the  pro¬ 
file  shape  as  derived  in  Eqs.  (2.31)  and  (2.32). 
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Therefore,  j6"  Is  the  only  quantity  which  sensitively  responds  to 
the  change  In  profile  shape  for  fixed  boundary  conditions.  In  order  to 
Investigate  how  related  geoacoustic  parameters  affect  this  quantity,  we 
shall  next  examine  the  analytic  expression  of  /J"  for  a  limiting  case, 
though  its  behavior  is  valid  for  a  large  portion  of  horizontal  ranges 
where  the  incident  angle  approaches  grazing. 


3.2.1  Geometric  Spreading  Near  Grazing  Incidence 
In  this  subsection  we  shall  demonstrate  from  the  profile  shapes  we 
have  chosen  that  d"  strongly  depends  on  the  sound  velocity  gradient  at 
the  top  of  the  layer  and  becomes  independent  of  the  profile  shape  near 
grazing  incidence. 

When  an  eigenray  of  a  primary  or  multiple  reflection  enters  the 
ocean  bottom  near  grazing,  it  turns  within  the  layer  at  relatively  shal¬ 
low  depth.  The  associated  phase  integral  is  calculated  by  integrating 
from  z  =  zt  to  z  =  0.  Let  us  consider  c^A'dp2  for  the  l/c2-l inear  pro¬ 
file  case 


t" 


2+  ZC 
2  3 
C0  y0 


+  2n 


d2n. 


d  p‘ 


z  =  0 


(3.19) 


where  the  first  term  corresponds  to  the  water  path  and  the  second  term 
corresponds  to  the  layer-path  given  by  Eq.  (3.4).  After  replacing  p  and 
Q  with  CjJzJ/cosaj  and  CjfzJ/sinaj  respectively,  we  rewrite  Eq.  (3.4)  as 


2 

=  ■  ajcprry 


2 

2COS  aj-l 
sinaj 


(3.20) 


where  is  the  grazing  angle  in  the  layer  at  z  =0.  We  next  di ffer- 
entiate  both  sides  of  Eq.  (3.1)  with  respect  to  z  to  obtain 


(2cos2arl).  (3.22) 

Since  we  are  considering  the  case  of  small  grazing  incidence,  the  small 
angle  approximations  to  the  leading  order  terms  may  replace  the  trigo¬ 
nometric  functions  and  Eq.  (3.22)  becomes 


(3.23) 

The  minus  sign  is  due  to  the  fact  that  positive  z  is  defined  to  be  up¬ 
ward  in  the  water  column,  and  therefore  the  velocity  gradient  in  the 
layer  is  negative. 

2  2 

Interestingly,  the  expressions  derived  for  d  n/dp  for  the 
Cj -linear  and  c^-1  inear  profiles  become  identical  to  the  one  given  by 
Eq.  (3.23)  when  small  angle  approximations  are  made.  Since  we  have 


fixed  cl  (0)  (and  c^  (h ) )  for  all  profile  shapes,  the  only  variable  is  the 


sound  velocity  gradient  at  the  ocean  bottom  for  the  rays  which 
completely  turn  within  the  layer. 

The  effect  of  this  Is  expected  to  be  more  critical  when  the  sum  of 
the  source  and  receiver  heists  is  sufficiently  small  so  that  the  first 
term  on  the  right  hand  side  of  Eq.  (3.19)  does  not  dominate  .  In 
order  to  quantitatively  determine  the  condition  where  the  water-path 
term  does  not  dominate  d",  we  shall  again  consider  the  case  of  small 
grazing  incidence  and  assume  that  Cq  s  c ^ (0 )  at  the  bottom.  Then, 
we  obtain  from  Eq.  (3.19)  that 

2  +  zQ  «  2n  a2c0/(dc1(0)/dz)  (3.24) 

where  a  is  the  grazing  angle  of  the  n-th  multiple  in  the  water  column. 

For  a  gradient  of  1  s'1,  the  right  hand  side  of  Eq.  (3.24) 
becomes  less  than  100  m  for  the  primary  reflection (n  =1).  As  the  order 
of  the  multiple  reflection  increases,  i.e.  n  in  Eq.  (3.19)  increases 
anq/or  the  sound  velocity  gradient  at  the  bottom  decreases,  the  layer- 
-path  contribution  increases.  Therefore,  a  larger  value  for  (z  +  zQ) 
satisfies  Eq.  (3.24). 

3.2.2  Geoacoustic  Models  with  Large  Layer  Thickness 

We  have  just  learned  that  geometric  spreading,  particularly  at 
long  ranges,  is  strongly  affected  by  the  sound  velocity  gradient  at  the 
bottom.  In  order  to  obtain  a  wide  range  in  the  sound  velocity  gradients 
among  the  three  profile  shapes,  while  keeping  the  same  boundary  condi¬ 
tions  at  both  interfaces,  one  must  have  a  relatively  large  layer  thick- 
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ness  and  sound  velocity  gradient. 

2 

For  the  given  endpoints,  the  sound  velocity  gradients  of  the  2/c^ 
ana  c^-linear  profiles  are  given  in  terms  of  the  endpoint- sound  ve¬ 
locities,  c1  (0 )  and  c^(H)  by 

p  dc,  (0)  cZA0)  ?  ? 

l/c‘-l  inear:  — ± -  — -  [c,  ( H)  -  c,  (0)]  (3.25) 

1  dz  2H  1  1 

■i  dc  i  ( 0 )  c  1  ( 0 )  -a 

cT-l inear:  -  -± -  Cc,J(H)  -  cf(0)3,  (3.26) 

1  dz  3H  1  1 

respectively ,  where  H  is  the  layer  thickness. 

We  next  define  a  "reference  gradient"  g  by  the  gradient  of  the 
c^-linear  profile  since  it  is  constant  throu^iout  the  layer.  Then, 

CjtH)  is  given  by 


c1(H)  =  Hg  ♦  0^(0).  (3.27) 

After  the  above  expression  for  Cj(H)  is  substituted,  Eqs.  (3.25)  and 
(3.26)  become 
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respectively.  Now  it  can  be  easily  shown  for  both  expressions  for  the 
gradient  d^  (0)/dz  that 


dcx(0) 

dz 

>  0 

as  g  >  0 

for  any  H, 

(3.30) 

dCj(O) 

t\y 

>  g 

as  H  >  0. 

(3.31) 

Therefore,  if  g  ano/or  H  become  small,  the  variations  in  the  velocity 
gradient  at  the  bottom  and  thus  the  entire  profile  shapes  become  close 
to  one  another. 

Fig.  3.1  shows  the  three  profile  shapes  between  fixed  end  points 

in  the  layer  whose  thickness  is  500  m  wi  th  g  =  1 .5  s'1.  For  this  ex- 

2  3 

ample  the  gradients  of  l/c£ ,  Cj  and  c^-1  inear  profiles  at  the  bottom 
are  0.8333  s'1,  1.5  s'1  and  2.375  s"1  respectively.  However,  if  the 
layer  thickness  is  reduced  to  50  m,  they  become  1.395  s'1,  1.5  s'1 
ana  1  .576  s'1,  respectively ,  for  the  same  reference  gradient  at  the 
bottom.  When  g  is  small,  the  differences  in  sound  velocity  gradient 
become  smaller  even  though  the  layer  thickness  is  increased  to  500  m 
(cf.  Fig.  3.2). 

One  interesting  non-dimensional  parameter  of  our  geoacoustic  model 
is  the  layer  thickness  divided  by  the  average  wavelength,  denoted  by 
H/<x1>.  A  large  layer  thickness  not  only  provides  substantial  varia¬ 
tions  among  profile  shapes  between  the  endpoints  but  also  sets  a  large 
value  for  this  parameter  for  reasonably  high  frequencies.  The  other 
extreme  case,  namely  the  layer  thickness  being  comparable  to  or  even 
smaller  than  the  wave! ength ,  is  also  of  interest  and  we  shall  examine 
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Fig.  3.1  Three  different  profile  shapes  for  g  =  1 .5  s-* 
l/c2,  (b )  c  and  (c )  c^-l  inear  types. 
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Fig.  3.2  Three  different  profile  shapes  for  g  =  0 
1/cS  (b )  c  and  (c  )  c^-l  inear  types. 
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geoacoustic  models  vrfiose  layer  thicknesses  are  varied  according  to  the 
size  of  this  quantity  in  Sec.  3.3. 

3.2.3  Geoacoustic  models 

While  we  shall  use  a  large  layer  in  order  to  study  the  effects  of 

different  profile  shapes,  we  shall  model  a  small  layer  with  a  single 

profile  type  since  the  exact  profile  shape  does  not  matter  much  for 

fixed  endpoint  veloci ties  and  a  small  layer  thickness.  We  shall  choose 
2 

die  l/c£-linear  profile  for  thin  layer  cases  since  for  this  profile 
the  Langer  reflection  coefficients  become  exact. 

Based  on  the  modeling  considerations  discussed  in  this  section,  we 
propose  studying  the  following  geoacoustic  models,  where  we  also  assume 
a  one  percent  sound  velocity  drop  at  the  top  of  the  layer,  a  phenomenon 
often  encountered  in  deep  ocean  bottom  experimental  datatl]: 

Model  I:  Thin  layer  model  with  a  5-m  layer  and  a  2  s"1  -gradient 

Model  II:  Small  layer  models  with  a  50 -m  layer  and 

a)  a  0.5  s'1 -gradient, 

b )  a  1 .5  s_1-gradient . 

Model  III:  Thick  layer  model  with  a  500-m  layer, 

a  .5  s"1- reference  gradient,  and 
2 

a)  a  1/c  -l  inear  profile, 

3 

b )  a  c  -1  inear  profile. 


Model  IV:  Thick  layer  model  with  a  500-m  layer, 
a  1 .5  s"*-reference  gradient,  and 
a)  a  1/c  -l  inear  profile 
b  )  a  c3-l  Inear  profll  e 

The  quantity  H/<x1>  is  0.1,  1  and  10  for  Models  I,  II,  and  III  and  IV, 
respectively . 

3.3  Reflected  Fields  for  The  Geoacoustic  Models 

We  still  have  to  fix  two  more  parameters  in  our  examples,  namely  the 
source  frequency  and  the  source-receiver  geometry.  Regarding  the  source 
frequency,  relatively  high  frequencies  must  be  used  in  order  for  the 
asymptotic  expansion  to  be  valid.  A  frequency  of  220  Hz  is  used  for  all 
of  the  examples  since  it  is  the  frequency  that  has  been  used  in  the 
fiela  experiments  by  Frisk[l]. 

The  value  of  the  sum  of  the  source  and  receiver  heights  off  the 
ocean  bottom  Z  is  defined  as 


This  has  also  been  varied  in  actual  field  experiments,  the  typical  range 
being  125-250  meters  from  the  ocean  bottom.  For  the  reflected  fields 
this  is  a  single  parameter.  If  the  direct  field  due  to  the  water  path 
is  Included  in  the  total  fields,  the  source  and  receiver  heights  must  be 
given  separately. 

The  reflected  fields  are  grouped  into  three  different  ray  types;  the 


! 

specular  reflection,  primary  reflection  and  multiple  reflections.  A 
schematic  of  these  rays  is  shown  in  Fig.  3.3.  The  total  reflected  field 
is  evaluated  from  the  coherent  sum  of  all  of  the  above  fields. 

In  order  to  obtain  the  reflected  fields  due  to  the  multiples  the 
infinite  series  of  Sommerfeld  integrals  for  the  multiples  may  be  trun¬ 
cated  at  the  m-th  term  so  that  the  following  condition  holds: 

m  CD 


2  RJ  s  2  RJ  =  rtotal;  0  <  ps. 

j  =0  j  =0  J 


(3.31) 


where  p$0  is  the  saddle  point  for  the  specularly  reflected  field  and 
psj  is  the  saddle  point  for  the  j-th  multiple.  This  is  a  more  rigor¬ 
ous  truncation  criterion  than  truncating  the  terms  whose  contributions 
are  below  a  certain  level  after  each  of  the  multiple  fields  is  computed 
because  it  guarantees  that  the  total  reflection  coefficient  in  the 
Sommerfeld  integral  is  well  approximated  before  the  integration  is  eval¬ 
uated. 


3.3.1  Model  I:  Geoacoustic  Model  with  Thin  Layer 
A  geoacoustic  model  with  a  5-meter  thick  layer  and  a  sound  veloci¬ 
ty  gradient  of  2  s"*(cf.  Fig.  3.4)  was  used  to  generate  the  reflected 
fields  for  Z  =  100  m  as  shown  in  Fig.  3.5.  Since  the  layer  thickness  is 
snail er  than  the  average  wavelength  in  the  layer (6. 77  m),  pure  ray  theo¬ 
ry  should  not  work  well.  Even  though  the  use  of  ray  theory  may  not  be 
valid  in  interpreting  the  reflected  fields,  the  ray  theory  based  expan¬ 
sion  of  the  total  reflection  coefficient  is  still  valid  and  we  can  still 


evaluate  each  of  the  Sommerfeld  integrals  for  the  f^'s. 


Specular 
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MODEL  I 


Fig.  3.4  A  thin  layer-model  with  a  5-m  layer. 
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Fig.  3.5  The  reflected  fields  for  the  Model  I  at  z 


Performing  this  analysis,  we  see  that  the  primary  reflection  domi¬ 
nates  the  total  reflected  field  up  to  a  horizontal  range  of  5000  m. 

Fig.  3.6  shows  a  ray  diagram  of  the  primary  reflections.  There  is  no 
caustic  generated  because  of  the  thin  layer.  The  specular  and  multiple 
reflections  are  of  about  equal  level,  though  both  are  substantially 
smaller  than  the  primary  reflection  for  most  ranges.  As  the  incident 
angle  approaches  grazing,  the  field  amplitudes  of  all  of  the  reflection 
types  become  almost  equal.  Since  the  differences  in  ray  path  between 
the  different  arrivals  are  very  small  for  the  chosen  source-receiver 
geomtry,  i.e.  the  water  column  is  twenty  times  thicker  than  the  bottom 
layer,  their  corresponding  saddle  points  are  very  close  to  each  other  at 
any  given  range.  This  is  confirmed  by  observing  that  the  specular  and 
multiple  returns  exhibit  minima  which  appear  at  the  same  incident  angle 
in  their  reflection  coefficient  plots  (similar  plots  are  shown  in  Figs. 
2.5  ana  2.6),  at  closely  located  ranges  near  800  m.  A  ray  diagram  of 
the  primary  reflections  is  shown  in  Fig.  3.6. 

We  conclude  from  the  above  observations  that  the  effect  of  the 
thin  layer,  even  with  a  rather  strong  sound  velocity  gradient  (such  as  2 
s'1  in  our  example),  is  weak  and  that  the  reflected  field  is  dominated 
by  the  subbottom  wi  th  which  the  primary  reflection  has  the  strongest 
interaction  except  at  ranges  associated  with  incidence  very  close  to 
grazing . 


3.3.2  Model  II:  Small  layer  thickness  with  variable  source- 
-  receiver  heists 

The  Models  II a  and  lib,  corresponding  to  two  different  sound  ve¬ 
locity  gradients,  0.5  s"*  and  1.5  s’*,  respectively (cf.  Fig.  3.7), 
are  considered  for  two  di fferent  val ues  of  source -receiver  height.  The 
reflected  fields  at  Z  =  50  m  and  Z  =  100  m  for  the  first  model  are  shown 
in  Fig.  3.8,  while  the  fields  for  the  second  model  for  these  values  of  Z 
are  also  shown  in  Fig.  3.9. 

For  both  models  highly  oscillatory  amplitude  variation  is  seen  at 
close  ranges  (r  <  50Qn).  This  is  due  to  the  strong  interaction  between 
the  specular  and  primary  reflections.  As  the  specular  reflection  ap¬ 
proaches  its  minimum  amplitude,  the  effect  of  the  multiples  becomes 
more  pronounced  and  alters  the  curve  of  the  total  reflected  field  which 
is  now  dominated  by  the  primary  reflection. 

Beyona  a  range  of  about  2000  m,  the  primary  reflection  decays  with 
range  because  the  amplitude  of  the  reflection  coefficient  for  the  prima¬ 
ry  reflection  starts  decaying  as  it  approaches  grazing  incidencefcf . 

Fig.  3 . 1Gb )  as  well  as  the  geometric  spreading  increases.  The  specular 
reflection  maintains  a  steady  level  due  to  a  sharp  increase  in  the  re¬ 
flection  coefficient  ampl i tude  (cf.  Fig.  3.10a)  countering  the  loss  due 
to  geometric  spreading.  In  both  models  the  primary  reflection  dominates 
for  most  ranges  at  Z  *  100  m. 

Hcwever,  at  Z  =  50  m  in  both  models,  the  specular  reflections 
become  dominant  over  the  primary  reflection  at  mid-ranges  (2500-3  500 
m).  This  is  due  to  the  fact  that  when  Z  is  reduced  to  50  m,  the  water- 
-path  contribution  to  the  geometric  spreading  becomes  less,  and  thus  the 
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Fig.  3.8b  The  reflected  fields  for  the  Model  1 1  a  at  z  +  z 
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3.9a  The  reflected  fields  for  the  Model  lib  at  z  +  z0  =  lOCm. 
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Fig.  3.9b  The  reflected  fields  for  the  Model  lib  at  z  * 
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contribution  of  the  layer-path  to  the  total  geometric  spreading  is  more 
emphasized. 

The  sharp  peak  in  the  primary  reflection  seen  at  range  1500  m  in 
Fig.  3.8b  (Model  Ila)  indicates  caustic  formation  near  the  ocean  bot¬ 
tom.  However,  as  Z  increases  to  100  m,  it  disappears.  This  disappear¬ 
ance  of  the  caustic  is  due  to  the  presence  of  the  subbottom  and  there 
exists  a  caustic  cut-off  between  these  depths.  Section  4.2  will  treat 
this  problem  in  more  detail  . 

Due  to  the  small  velocity  gradient  in  Model  Ila  the  field  contrib¬ 
ution  from  the  multiples  is  small  even  at  long  ranges  (r  >  3000  m). 
However,  when  the  velocity  gradient  is  increased  to  1.5  s-1,  as  shown 
in  Model  lib,  their  contribution  becomes  one  of  the  major  ones  even  at 
relatively  close  ranges  for  small  Z  (cf.  Fig.  3.9b).  Therefore,  a  large 
sound  velocity  gradient  at  the  bottom  and  a  small  Z  contribute  to  a 
larger  multiple  fiel d  ampl i tude . 

3.3.3  Model  III:  Thick  Layer  with  Small  Gradient 

2 

Two  different  profile  shapes  are  consider  in  this  model:  1/c 

3 

ano  c  -linear(cf.  Fig.  3.11).  As  discussed  in  Sec.  3.2.2  and  shown  in 
Fig.  3.2,  a  small  reference  gradient(.5  s"*in  this  model)  creates  only 
snail  differences  in  the  profile  shapes.  However,  even  a  small  differ¬ 
ence  in  profile  shape  can  make  large  differences  in  the  geometric 
spreading  at  long  ranges.  The  large  layer  thickness  makes  it  possible 
to  forma  caustic  extending  into  the  water  column(cf.  Fig.  3.12).  In 
this  case  the  total  reflected  fields  are  characterized  by  the  fields 
beyond  the  caustics  formed  by  the  primary  and  mul  ti  pi  es  (cf .  Fig.  3.13), 


•:V/ 


o  O  □  SPECULAR  REFLECTION 
0-0  O  PRIMARY  REFLECTION 
A  A  A  MULTIPLE  REFLECTIONS 
+-+-+  TOTAL  REFLECTION 

-  1  -  I  ->  ■  i  I — »  >  -I  -I  I  »  "  i  -  I — i  [  i  1  i  — >  — ;■  "I r  - t  — i  — I 

0  1000  2000  3000  4000  5000 

RANGE  (M) 


The  reflected  fields  for  the  Model  Ilia  at  z  +  zq 


L_J 

3 

t; 

2  ‘10° ' 
< 

UJ 

> 

h-  -120- 
< 

UJ 

I  1 

o  c 

£  <j 

y 

B  -O . O  SPECULAR  REFLECTION 

o— e~e  PRIMARY  REFLECTION 

a-a-a  multiple  reflections 
+-+-+  TOTAL  reflection 

-  IUU  \ 

c 

I  1000  2000  3000  4000  5000 

RANGE  (M) 

Fig.  3.13b 

The  reflected  fields  for  the  Model  1 1 1 b  at  z  ♦  zq 

IIII1 


[I] 


while  the  overall  field  amplitudes  of  all  reflection  types  are  very 
similar  otherwise. 

These  strongly  interfering  primary  and  multiple  reflection  ampli¬ 
tudes,  which  result  from  caustic  formations,  are  due  to  multi-saddle 
point  contributions (cf.  Sec.  2.3.2).  The  difference  in  locations  of  the 
appearance  of  these  features  are  very  clear  in  both  profile  types. 
However,  the  patterns  of  the  amplitudes  of  both  primary  and  multiples 
are  almost  identical  in  both  profile  types.  The  similarity  in  the  field 
amplitudes  near  the  caustic  will  be  discussed  in  Sec.  4.3. 

3.3.4  Model  IV:  Thick  Layer  with  Large  Gradient 

The  sound  velocity  gradient  at  the  top  of  the  layer  varies  strong¬ 
ly  in  this  model,  given  that  the  boundary  conditions  at  both  of  the 
interfaces  are  the  same(cf.  Fig.  3.14).  As  shown  in  Figs.  3.15  and 
3.16,  the  caustic  formations  of  the  primary  reflections  are  quite  dif¬ 
ferent  in  the  1/c2  and  c3-linear  layer  cases.  In  each  case  the 
interference  patterns  due  to  multi-saddle  points  characterize  the  be¬ 
haviors  of  the  primary  and  multiple  refl ections(cf .  Fig.  3.17).  In  the 
l/c^-linear  case  shewn  in  Fig.  3.15  the  caustic  has  two  branches  and  a 
cusp  at  the  point  where  they  meet.  The  interference  of  the  primary 
reflections  only  takes  place  within  the  ranges  surrounded  by  the  caustic 

3 

brandies.  In  the  case  of  the  c  -linear  profile,  there  is  only  one 
caustic  branch  which  extends  into  the  water  column  up  to  the  caustic 
cut-off  point  (cf.  Sec.  4.2).  For  large  layer  thickness,  the  caustic 
branch  can  extend  to  large  distances  above  the  ocean  bottom.  Thus  the 
interference  of  the  primary  reflections  can  extend  to  longer  horizontal 
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Fig.  3.14  A  thick  layer  model  wi th  a  500-m  layer  where  g  =  1 .5  S“*. 

(a)  The  1/ c 2- 1  Inear  profile  and  (b)  the  c^-l  inear  profile. 
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Fiy.  3.15  A  ray  diagram  of  the  primary  reflections  for  the  Model  IVa 
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ranges  until  the  subbottom  reflection  terminates  this  behavior (cf.  Fig. 
3.16). 

Comparing  the  reflected  fields  generated  from  both  cases  with  dif¬ 
ferent  sound  velocity  gradients,  we  observe  that  the  total  field  from 
Model  IVb(cf.  Fig.  3.17b),  which  has  a  larger  gradient  at  the  bottom 
than  Model  IVa(cf.  Fig.  3.17a),  shows  lower  amplitude  levels  beyond  the 
primary  caustic  range.  One  may  expect  from  our  discussion  in  Sec.  3.2.1 
that  the  primary  reflection  for  a  geoacoustic  model  with  a  larger  sound 
velocity  gradient  should  have  a  larger  amplitude.  However,  for  our 
typical  source-receiver  geometries,  i.e.  Z  >  100  m,  the  condition  given 
by  Eq.  (3.24)  is  violated  and  the  water-path  contribution  to  the  total 
geometric  spreading  becomes  substantially  larger  than  that  of  the  layer 
as  the  eigenray  approaches  grazing  incidence. 

Since  OgCcf.  Eq.  (3.24)]  of  the  primary  eigenrays  reaching  long 
ranges  in  Model  IVb  are  smaller,  an  even  greater  water-path  contribution 
is  added  to  the  total  geometric  spreading.  As  a  result,  the  primary 
reflection  from  Model  IVb  has  less  amplitude  than  that  of  Model  IVa. 
Furthermore,  since  the  primary  reflections  are  the  dominant  contributor 
an o  control  the  overall  level  of  the  total  reflected  fields  for  both 
cases,  the  total  reflected  field  from  Model  IVb  shows  less  relative 
amplitude  even  with  the  greater  velocity  gradient  near  the  bottom. 

Multiples,  on  the  other  hand,  effectively  reduce  geometric  spread¬ 
ing  due  to  both  the  water-path  and  the  layer-path,  thus  multiple  con¬ 
tribution  becomes  larger  from  the  geoacoustic  model  with  a  larger  sound 
velocity  gradient  at  the  bottom.  The  multiple  reflections  also  form 
similar  Interference  patterns  to  the  primary  reflections  and  are  respon- 
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Fig.  3.17a  The  reflected  fields  for  the  Model  IVa  at  z  ♦  zq  =  100 
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s ibl e  for  the  highly  oscillatory  patterns  in  the  total  reflected  fields 
observed  at  long  ranges  beyond  the  primary  caustic  ranges  for  both  mod¬ 
els. 
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Chapter  4:  Analytic  Studies  on  Caustics 

Due  to  the  limitations  of  the  method  of  steepest  descents,  the 
analytic  evaluation  of  the  reflected  fields  near  caustics  was  excluded 
in  Chapter  2.  However,  understanding  the  parametric  dependence  of  the 
caustic  shapes  and  the  strength  of  the  field  near  the  caustics  can  help 
recover  geoacoustic  properties  of  the  ocean  bottom[2].  In  this  chapter 
we  shall  focus  on  this  subject  and  discuss  the  caustic  behavior,  pri¬ 
marily  in  a  qualitative  manner,  for  the  proposed  sound  velocity  profile 
shapes . 

4.1  The  caustic  equations  for  the  proposed  geoacoustic  models 

Brekhovskikhtl]*  derived  the  equations  of  caustics  formed  in  a 
homogeneous  upper  half-space  overlying  an  Inhomogeneous  halfspace  for 
the  n2-linear  and  c-linear  profiles  where  the  sound  velocity  is  con¬ 
tinuous  across  the  boundary.  This  section  is  an  extension  of  his  work 
to  derive  a  caustic  equation  for  the  case  of  the  c  -linear  profile 
which  provides  a  more  realistic  profile  shape.  The  equations  of  cau$^ 
tics  for  all  three  cases  are  derived  for  the  n-th  multiples  and  there¬ 
fore  they  are  more  general  forms  than  those  of  Brekhovsk ikh .  In  this 
section  we  shall  neglect  the  effect  of  the  subbottom  on  caustic  formation 

*  The  results  given  in  Ref.  1,  Eqs.  (46,21)  through  (46.23)  appear 
to  have  errors.  Our  corrected  formulas  are  given  in  Sec.  4.1.2 
of  this  chapter . 


and  will  reserve  a  discussion  of  that  subject  for  Sec.  4.2. 


The  cases  treated  In  this  chapter  are  limited  to  the  geoacoustic 
models  where  there  Is  no  discontinuity  in  sound  velocity  at  the  ocean 
bottom  interface.  Even  though  we  have  chosen  less  comprehensive  geo¬ 
acoustic  models  than  those  we  have  used  earlier  in  order  to  make  the 
algebra  easier,  we  may  still  learn  a  great  deal  about  the  caustics  from 
these  special  cases  for  the  various  inhomogeneities  in  the  layer. 


4.1.1  The  caustic  equations  for  the  c  -1  inear  profile 
The  conaitions  for  the  caustic  in  the  water  column  are  given  by 
Eqs.  (2.40)  and  (2.44)  when  p$  is  such  that  the  left  hand  side  of  Eq. 
(2.44)  vanishes.  After  Eqs.  (3.15)  and  (3.16)  are  evaluated  at  the 
water-bottom  interface  and  substituted  in  Eqs.  (2.40)  and  (2.44)  respec¬ 
tively,  we  obtain 


—  (z  ♦  zn)  ♦  r  ♦  (i  .  p?b|/3)1/2(2p"2  ♦  b|/3)  x  0  (4.1) 

r0  u  c  a3pc  c  J  c  j 


z  *  z 


0  6n  , ,  2^2/3%, .  2. 2/3.-1/2  n 
zrr  ■  7“T (2  •  p?3  )(1  •  pcb3  }  =0 

Vo  a3pc 


(4.2) 


where  pc  is  the  horizontal  slowness  which  satisfies  the  above  equa¬ 
tions,  a3  ana  are  the  coefficients  in  Eq.  (3.13),  and  r£  is  the 
caustic  range.  If  the  sound  velocity  in  the  lower  half-space  is  given  by 


c3 (z  >  »  a3z  ♦  cj 


(4.3) 
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.*  /  -■ 


K  : 


where  cQ  is  the  sound  velocity  In  the  upper  half-space,  Eqs.  (4.1)  and 
(4.2)  may  be  reduced  to 


rc  *  £  <2  *  20>  -  ifj£  11  -  £o)1/2,2Pc2  *  c0>  I4-4' 


0  ^  6n  /0  _2„2W,  _2_2x-l/2  „ 

___  ♦  — T  (2  -  pcc0)(l  -  pcc0)  =  0 

VO  a3PC 


(4.5) 


Rewriting  yq  an d  pc  in  terms  of  the  grazing  angle  a  in  the  me¬ 
dium  containing  the  inci dent  wave,  we  then  obtain 


2c]p  2 

r.  =  (z  +  Zq)co ta  -  tana  ( — “  +  U 

3  COS  a 


(4.6) 


7+_  6nc0  sinZa  (9  2. 

z+  2o*  "  TT  (  * cos  a) 

3  COS  a 

Solving  for  cos2a  from  Eq.  (4.7),  we  obtain 


(4.7) 


COS2a 


3nC0  M? 

- « - *L3M1 - ^(z^zn))1/2]  (4.8) 

a3(z  +  20}  +  6nc0  3nc0 


Since  a^  <  0,  there  is  no  limit  on  z  +  Zq  and  the  caustic  branch  extends 
to  infinity  in  the  upper  halfspace.  Since  cos2a  <  1 ,  we  must  choose 
the  minus  sign  for  the  square  root  for  any  a^(z  +  Zq). 

Therefore,  the  location  of  the  caustic  may  be  obtained  for  a  given 
z  ana  Zq  by  substituting 
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a  *  Arccos 


(4.9) 


t 

I 

i 

i 


3nc: 


a3(z  *  z0)  *  6ncj 


[3  -  (1  - 


4  a 


3nc 


3  (z  .  z,))1'2]1'2 


In  to  Eq.  (4.6)  In  order  to  obtain  the  caustic  range  r£.  Fig.  4.1a  shows 

a  plot  of  rc  vs.  z  for  a  source  hei<Jit  of  50  m,  Cq  =  1500  m/s,  n  =  1 

6  2  3 

(primary  reflection)  and  a3  =  -6.75  x  10  m /s  . 


4.1.2  The  caustic  equation  for  the  1/c^-l  inear  profile 

2 

The  caustic  equation  in  the  water  column  for  the  1/c  -linear  pro¬ 
file  may  be  similarly  derived  and  it  is  expressed  by  the  following 
equations: 


f 

i  ,* 

K 

i 

■ 


a  =  Arcs 


Injjy  [1.(1-  2>2<Z  *  ;°^)1/2]1/2 


4n  _.._2 

7 


r_  =  — sin  a  cos  a 


(4.10a) 

(4.10b) 


Vo 


where  a  is  the  grazing  angle  of  the  ray  in  the  water  column  and  a^  is 

the  coefficient  of  z  in  Eq.  (3.1).  For  a  given  z  and  zQ,  one  may 

obtain  a  from  Eq.  (4.10a ) ,  and  substi  tute  i  t  in  Eq.  ( 4 .10b  )  to  obtain 

the  caustic  range  rQ .  Since  a  has  two  values  for  a  given  z,  zQ  and 

aj ,  then  rQ  does  also.  This  is  a  very  important  characteristic  of 

2 

the  caustic  shape  generated  by  the  1/c  -linear  profile  whose  caustic 
branches  never  extend  to  positive  infinity  in  the  upper  halfspace. 
Instead  they  form  a  cusped  caustic  known  as  a  butterfly  caustic (cf. 
Figs.  4.  2a,b ). 
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At  the  cusp  a  has  a  single  value  satisfying  the  condition. 


2ax(z  +  zQ)  c£/n  =  1 
Thus,  we  obtain  from  Eq.  (4.10) 


(4.11a) 


a  =  ir/6  , 


rc  ■ 


33/2n 

2^ 


(4.11b) 


Equation  (4.11)  indicates  that  the  grazing  angle  for  the  cusp  is  inde¬ 
pendent  of  z  +  Zq  and  aj  ,  and  is  always  30  degrees.  The  regions  of 
a  which  correcpond  to  separate  branches  of  the  caustic  in  the  water 
column  are  given  as 


on  the  branch  2-3:  aj  <  a  <  it/6 
on  the  branch  1-2:  ir/6  <  a  < 

where  ^  -  Arcsin  -^[1  ^  (1  -  Za^ZgCg/n)  ^  ]  (4.12) 

4.1.3  The  caustic  equation  for  the  c-linear  profile 
The  derivation  of  the  caustic  equation  for  the  c-linear  profile  is 
straightforward  and  the  result  is  given  by 


z  ♦ 


(4.13) 
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lere  a2  Is  the  coefficient  of  z  in  Eq.  (3.7), 


4.1.4  Caustic  equations  for  small  grazing  angl es 
Fig.  4.3  shows  the  curves  of  the  caustics  in  the  water  column  for 
the  three  different  profile  shapes  given  that  all  of  their  sound  veloc¬ 
ity  gradients  at  the  ocean  bottom  interface  are  fixed  at  1  sec-1  (cf. 
Fig.  4.4),  Notice  that  all  of  them  have  very  similar  curves  (consider- 

O 

ing  tlie  upper  branch  of  the  1/c  -linear  profile)  especially  near  the 
ocean  bottom.  This  is  due  to  the  fact  that  all  of  the  caustic  equations 
become  identical  as  the  grazing  angles  become  very  small.  If  small 
angle  approximations  are  substituted  in  the  pairs  of  caustic  equations, 
Eqs.  (4.6)  and  (4.7),  and  Eqs.  (4.10a  ,b),  they  can  be  reduced  to  (Eq. 
(4.13)  is  exact): 


p 

for  the  1/c  -linear  profile: 


a  C 

z  ♦  zQ  s  /  (upper  branch)  (4.14) 


for  the  c  -1  inear  profile: 


z  +  z  s  - 


(4.15) 


where  a  is  approximated  by 


,  .  r(Z  +  2<J)C>2 

*  *  £'“2*o 


(4.16) 


dcf  (0) 


i  *  1,3 


For  the  same  gradient  at  the  ocean  bottom  (z  =0),  a^  and  a^  may  be 
given  in  terms  of  a2  as 


If  a^  ana  a^  are  substituted  in  Eqs.  (4.14)  and  (4.15)  respec¬ 
tively,  they  may  be  identically  reduced  to  Eq.  (4.13)  for  small  grazing 
angles  associated  wi  th  caustics  formed  near  the  ocean  bottom.  Since 
small  a  corresponds  to  small  c'Zg  [cf.  Eq.  (4.16)],  we  can  conclude  that 
the  caustic  curves  near  the  ocean  bottom(z  is  also  small)  are  controlled 
by  the  gradient  of  the  profile  near  the  bottom,  not  the  overall  shape  of 
the  profile.  This  observation  is  particularly  important  when  the  sub¬ 
bottom  in  present  since  the  caustics  in  the  water  column  are  formed  only 
near  the  bottom  for  this  case,  which  is  discussed  in  the  next  section. 

4.2  The  effects  of  the  subbottom  interface  on  the  caustic 

in  the  water  column 

In  Section  4.1  the  inhomogeneous  medium  was  treated  as  an  infinite 
halfspace  for  which  caustics  are  formed  by  purely  refracted  rays.  How¬ 
ever,  the  presence  of  a  subbottom  reflector  in  the  geoacoustic  models 
under  consideration  serves  to  eliminate  rays  with  certain  ranges  of  ho¬ 
rizontal  slowness  p  from  caustic  formation.  Brekhovskikh  has  delineated 
the  regions  of  slowness  for  different  branches  of  the  caustics  for  the 
n  -linear  and  c-linear  cases.  The  presence  of  the  subbottom  simply 


sets  the  limits  on  the  regions  of  slowness  forming  the  various  caustic 
branches,  but  does  not  create  additional  branches. 

O 

For  the  case  of  the  1/c  -linear  profile,  the  caustic  branches  will 
be  eliminated  in  numerical  order  as  labeled  in  Fig.  4.2a  as  the  rays 
forming  the  caustic  start  reflecting  off  the  subbottom  interface.  For 
the  other  two  profile  shapes,  which  share  similar  caustic  curves  with 
only  two  branches,  the  caustic  branch  extending  into  the  water  column 
away  from  the  bottom  will  be  eliminated.  The  order  of  disappearance  of 
the  caustic  branch  is  again  indicated  by  the  numbers  labeled  along  the 
caustic  branches  in  Fig.  4.1a. 

Thus,  the  branch  of  the  caustic  extending  into  the  water  column  is 

cut  off  as  the  last  caustic -forming  refracting  ray  interacts  with  the 

subbottom.  The  location  of  the  caustic  cut-off  may  be  easily  obtained 

by  substituting  the  p  of  the  last  turning  ray  into  the  ray  equation,  say 
3 

Eq.  (4.1)  for  c  -linear  profile,  to  obtain  r  ,  and  then  substituting 

V 

r  into  the  equation  of  the  caustic  to  find  the  receiver  depth  z  for  a 
given  zQ.  Figs.  4.5a,b  show  ray  trace  diagrams  with  and  without  the 
subbottom  for  the  c-linear  profile  in  the  lower  half-space.  Fig.  4.5a 
clearly  shows  the  caustic  cut-off  in  the  water  column  where  the  rays 
start  radiating  into  a  wedge-shaped  region  without  crossing  each  other. 

4.3  The  fields  near  caustics 

As  discussed  in  Chapter  2,  the  asymptotic  solution  derived  in  that 
chapter  is  not  valid  for  the  fields  near  the  caustics  vrfiere  0"  in  Eqs. 
(2.41)  ana  (2.42)  becomes  very  small  and  vanishes  completely  at  the 
caustics.  The  fields  in  these  regions  can  be  evaluated  by  numerical 
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4.5  Ray  diagrams  of  the  primary  reflection 
profile  (a)  with  the  subbottom  interfa 
subbotm  interface. 
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Integration  methods.  Analytic  solutions  to  such  fields  also  exist  with 
some  limitations.  Brekhovskikh 's  caustic  correction  to  ray  theory  [3] 
is  only  valid  for  the  fields  very  near  the  caustics.  A  uniform  asymp¬ 
totic  solution  such  as  the  one  by  ifiite  and  Pedersen[4]  continuously 
valid  across  the  caustic  involves  a  great  deal  of  mathematics,  and  it  is 
rather  difficult  to  obtain  any  physical  interpretation  of  the  fields. 

In  this  section  we  shall  emphasize  the  physical  interpretation  of 
the  various  caustics  resulting  from  the  different  profile  types  discuss¬ 
ed  in  Sec.  3.1.  The  shapes  and  locations  of  caustics  in  the  water  col¬ 
umn  have  been  discussed  in  Secs.  4.1  and  4.2.  Our  main  interests  in 
this  section  are  to  evaluate  analytically  the  strength  of  the  caustic 
ana  geometric  decay  of  the  fields  near  the  caustic  as  one  moves  away 
from  the  caustic  with  horizontal  range.  Therefore,  analytic  solutions 
derived  by  Brekhovskikh  which  are  only  valid  near  the  caustic  serve  our 
neeas,  ana  we  shall  discuss  the  above  problems  in  terms  of  the  three 
profile  types  having  the  same  sound  velocity  gradients  at  the  bottom(cf. 
Fig.  4.4). 

We  have  modified  Brekhovskikh 's  analytic  solution  to  include  the 
Langer  reflection  coefficient  R n (p )  to  the  fields  near  the  first  order 
caustic  (i  .e.  dtf/dp  -  d2b/dp2  =  0  but  d3j6/dp3  /  0  at  p  =  pc) 
to  obtain 


5/6  1/6  iir/4,pc»  1/2  Rn(pc}  U3 


2J/  V'ue  ,ir/^(-^) 
r 


Ai  (q) 


(4.17) 
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where  Ai  is  the  Airy  function  whose  variable  q  is  given  by 


1/3  A<PC>  _1/3 

q  =  s9n[/>,,,{pc)]  21/Ju  - (r  -  r  )  (4.18) 

dpJ 

The  value  of  pc  varies  depending  on  the  profile  shapes  and  the 
source-receiver  or ientations .  Thus,  the  coefficient  of  the  Airy  func¬ 
tion  in  Eq.  (4.17)  involves  many  parameters  to  be  examined.  However,  as 
discussed  in  Sec.  4.1.4,  if  the  caustic  fields  of  interest  are  relative¬ 
ly  near  the  bottom,  the  variations  in  pc,  r  and  1^  are  very  small  since 
the  profile  shapes  at  the  top  of  the  layer  are  very  similar.  Therefore, 
the  most  sensitive  quantity  is  which  also  plays  an  important 

role  as  the  weighting  factor  in  the  argument  of  the  Airy  function,  which 
dominates  the  variation  with  range  of  the  field  amplitude  near  the  caus¬ 
tic. 

In  order  to  evaluate  the  relative  strength  of  the  caustic  and  the 
decay  with  horizontal  range,  it  is  useful  to  overlay  the  Airy  function 
multi plied  by |d3/6/dp3('1^3  evaluated  at  the  pc  corresponding  to  each 
profile  shape,  and  plotted  against  the  horizontal  distance  away  from  the 
caustic.  For  two  different  values  of  the  source  and  receiver  heights,  z 
♦  Zg  s  250  m  and  50  m,  Figs.  4.6a  ,b  show  two  sets  of  curves  of  the  men¬ 
tioned  quantity  for  the  three  different  profile  types  based  on  the  geo¬ 
acoustic  models  shown  in  Fig.  4.3.  The  frequency  of  the  harmonic  point 
source  is  220  Hz. 

For  Zq  =  50  m,  at  z  =  0  m  all  three  caustic  curves  are  very 
similar  (cf.  Fig.  4.3)  and  the  quanti  ties  of  I  b"'| '^3Ai  in  Fig.  4.6b 
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also  show  close  simil arli ties  among  the  three  profile  types.  On  the 
other  hand,  at  z  =  200  m  |<J'"|'1/^3Ai  for  the  c  and  c3-l inear  profiles 
similarly  have  greater  amplitudes  and  spatial  oscillations  than  at  z  =  0 

p 

m,  but  the  behavior  of  the  1/c  -linear  profile  remains  almost  the  same 
at  both  of  the  depths.  Even  though  the  caustic  curves  in  Fig.  4.3  show 
well  separated  traces  near  z  =  200  m,  one  can  observe  a  similarity  bet- 
ween  the  c  ana  c  -linear  profiles  in  terms  of  the  radius  of  curvature 
while  the  curve  for  the  l/t^-linear  profile  yields  a  sharper  turning. 

The  difference  in  the  field  strength  can  be  physically  explained 
by  ray  interpretations  as  follows.  Let  us  consider  rays  forming  two 
types  of  caustics  in  a  homogeneous  medium.  One  of  them  is  formed  along 
a  curve  of  a  small  radius  of  the  curvature  and  the  other  is  formed  along 
a  large  radius  of  the  curvature  as  shown  in  Figs.  4.7a  and  4.7b,  respec¬ 
tively.  As  seen  in  both  figures,  the  rays  forming  the  caustic  with  the 
larger  radius  of  the  curvature  are  more  densely  interacting  wi th  each 
other  near  the  caustic.  This  results  in  larger  and  highly  oscillatory 
interference  field  amplitude  patterns  near  the  caustic  as  observed  in 
Fig.  4.6a  for  the  c  and  c  -1  inear  profiles. 

We  can  conclude  from  the  above  observations  that  when  the  caustic 
curves  are  similar,  the  associated  fields  near  such  caustics  are  also 
similar  and  they  are  more  directly  affected  by  the  curvatures  of  the 
caustic  curves  than  by  the  shapes  of  the  profiles. 
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Chapters:  Analysis  of  Experimental  Data 

We  have  examined  various  geoacoustic  models  using  an  asymptotic 
expansion  method  based  on  the  Langer  reflection  coefficients  in  Chapter 
3.  In  this  chapter  we  shall  examine  two  sets  of  experimental  data  ob¬ 
tained  from  field  experiments  and  estimate  the  geoacoustic  properties  of 
the  ocean  bottom  from  these  data  using  the  asymptotic  expansion  method. 
After  finding  geoacoustic  models  such  that  the  resulting  total  fields 
agree  with  the  measurements  satisfactorily,  we  shall  demonstrate  the 
energy  partitioning  of  the  eigenrays  with  horizontal  range,  so  that  the 
physical  mechanisms  causing  the  real  fields  may  be  understood.  We  shall 
also  discuss  the  effect  of  source/receiver  geometry  on  the  determination 
of  geoacoustic  models. 

Data  Setl:  Icelandic  Basin 

Field  measurements  were  taken  in  the  Icelandic  Basin  (59  31.0'N, 

17  20.8 'W)  by  Frisk  and  his  co-workers[l  J.  The  total  acoustic  pressure 
field  due  to  a  220-Hz  harmonic  point  source  towed  at  an  elevation  of  157 
m  above  the  ocean  bottom(cf.  Fig.  1.2)  was  measured  by  a  receiver  moored 
at  54.6  m  above  the  bottom.  It  has  been  found  by  FriskCl ](with 
modifications  by  the  author)  that  the  geoacoustic  model  shown  in  Fig. 

5.1  seems  to  provide  the  best  synthetic  fields  which  match  with  these 
measurements.  Both  synthetic  and  real  data  are  shown  in  Fig.  5.2. 

We  next  generate  reflected  fields  based  on  this  geoacoustic  model 
and  group  them  into  three  ray  types:  specular,  primary  and  multiple 
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A  comparison  between  experiment  and  asymptotic  method  using 
the  geoacoustic  model  shown  in  Fig.  5.1.  The  source 
frequency  is  220  Hz  and  the  source  and  receiver  heights  are 
157  m  and  54.6  m,  respectively. 


reflections.  As  shown  in  Fig.  5.3,  the  primary  reflection  clearly  domi¬ 
nates  tli e  total  reflected  field  (excluding  the  direct  field)  for  ranges 
up  to  5  km  even  though  the  sound  velocity  gradient  is  the  largest  we 
have  examined  so  far  and  one  would  expect  a  greater  contribution  from 
the  mul  ti pi  es (cf .  Sec.  3.2.1).  This  is  due  to  the  fact  that  the  water 
column  occupies  a  large  portion  of  the  total  ray  path  even  for  the  mul- 
tipleslcf.  Fig.  5.4),  and  the  condition  given  by  Eq.  (3.24)  is  violated 
(the  R.  H.  S.  of  Eq.  (3.24)  is  about  30  m,  while  (z  +  zQ)  is  211.6  m). 

Therefore,  although  the  geometric  spreading  due  to  the  layer-path 
is  relatively  small  for  the  steep  sound  velocity  gradient  at  the  bottom, 
this  does  not  contribute  much  to  reducing  the  total  geometric  spreading 
of  the  multiples.  Given  an  insignificant  contribution  to  the  field 
amplitude  from  the  steep  sound  velocity  gradient  at  the  bottom,  the 
primary  reflection  which  has  the  maximum  reflection  coefficient  ampli¬ 
tude  after  critical  incidence  dominates  the  total  reflected  field. 

Based  on  these  observations  we  can  also  conclude  that  the  sound 
velocity  at  the  top  of  the  layer  and  the  gradient  in  the  layer,  which 
most  affect  the  behavior  of  the  primary  reflection,  particularly  its 
geometric  spreading,  are  mainly  responsible  for  the  nature  of  the  total 
reflected  field  at  these  ranges.  The  geoacoustic  model  shown  in  Fig. 

5.1  is  the  result  of  using  such  modeling  criteria  to  determine  the  sound 
velocity  at  the  top  of  the  layer  and  the  gradient  in  the  layer  to  im¬ 
prove  the  agreement  between  the  measurement  and  synthetic  fields. 

The  fact  that  in  this  case  geometric  spreading  is  heavily  depen¬ 
dent  upon  the  water-path  and  not  the  layer-path  makes  the  identification 
of  geoacoustic  models  somewhat  more  difficult,  particularly  the  shape  of 
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The  reflected  fields  for  the  geoacoustic  model  shown  in  Fig 
5.1 .  at  z  +  Zg  s  211 .6  m 


Fig.  5.4  Ray  diagrams  of  the  (a)  primary  reflections  and  (b )  second 
multiple  reflections  for  the  model  shown  in  Fig.  5.1. 


.  5.5  Synthetic  data  for  the  (a)  l/c2-l  inear  profile  and  (b) 

c^-linear  profile  instead  of  the  c  linear  profile  in  the 

layer  shown  in  Fig.  5.1.  • 
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the  profile  associated  with  the  sound  velocity  gradient  at  the  bottom. 
This  is  n£t  the  same  as  saying  that  matching  the  synthetic  fields  to  the 
measurements  is  difficult. 

As  an  example.  Fig.  5.5  shows  two  sets  of  the  total  field:  one 

2  ^ 
from  the  1/c  -linear  profile  and  the  other  from  the  c  -linear  pro¬ 
file  for  fixed  endpoints  of  the  sound  velocity  in  the  layer.  It  is 
evident  that  all  of  the  total  fields  generated  from  the  geoacoustic 
models  with  different  profile  types,  including  the  one  shown  in  Fig. 

5.2,  have  good  agreement  with  the  measurements. 

Therefore,  when  the  source  and  receiver  are  placed  at  relatively 
high  elevation  above  the  bottom  so  that  the  water-path  dominates  the 
total  geometric  spreading,  the  total  fields  measured  at  these  receiver 
locations  are  insensitive  to  the  shape  of  the  bottom  layer  profile,  and 
thus  the  identification  of  this  property  of  the  bottom  structure  becomes 
very  difficult.  In  the  next  example,  we  shall  demonstrate  that  the 
effect  of  the  profile  shape  is  indeed  significant  when  the  source  and 
receiver  heights  are  reduced. 

Data  Set  2:  Hatteras  Abyssal  Plain 

We  next  examine  experimental  data  collected  in  the  Hatteras 
Abyssal  Plain  (34  00.0' N,  67  00.0'W)  with  the  same  experimental 
configuration  as  in  the  previ-  ous  example  except  that  the  source  and 
receiver  heights  are  now  180  and  2.7  meters,  respectively.  The 
estimated  geoacoustic  model  for  this  region  is  shown  in  Fig.  5.6.  where 
the  average  water  sound  velocity  was  used  from  the  model  by  Frisk[2]. 

The  measurement  of  the  total  field  is  shown  along  with  the  synthetic 
field  based  on  our  geoacoustic  model  in 
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Fig.  5.7.  Due  to  the  comparable  size  of  the  source  plus  receiver 
heights  and  the  sediment  layer  thickness{c ,f .  Ray  trace  diagram  in  Fig. 
5.8),  this  example  demonstrates  a  sensitive  dependence  on  the  shape  of 
the  sound  velocity  profile  in  the  layer  despite  the  small  reference 
sound  velocity  gradient (dc/dz  =  1  s'*  as  opposed  to  a  value  of  nearly 
3  s"1  for  the  previous  example). 

Except  at  ranges  near  r  =  2000  m,  both  real  and  synthetic  fields 
agree  well.  This  agreement  is  especially  good  at  very  long  ranges  as 

3 

the  result  of  the  c  -linear  profile  being  used  for  the  sound  velocity 

in  the  layer.  If  one  instead  takes  the  profile  in  the  layer  to  be 
2 

1/c  -linear  for  the  same  fixed  endpoints,  the  resulting  synthetic 
field  shows  a  very  di fferent  behavior ,  as  seen  in  Fig.  5.9.  We  note 
that  the  large  peak  at  and  just  beyond  the  caustic  has  become  narrower. 
As  shown  in  Fig.  5.10,  this  is  due  to  there  being  a  narrower  horizontal 
range  band  in  which  the  primary  reflections  can  interfere  with  each 
other. 

The  profile-dependent  differences  in  the  reflected  field  ampli¬ 
tudes  seen  at  long  range  can  be  analyzed  by  comparing  their  grouped 
reflected  fields,  shown  in  Figs.  5.11  and  5.12  for  the  C*  and  1/c2- 
-linear  profile  shapes,  respectively.  Except  near  the  caustic  ranges, 
the  gross  variations  of  the  field  amplitudes  of  each  reflection  are 
almost  the  same.  However,  the  total  reflected  field  amplitude  behaves 
somewhat  differently  in  each  case  with  respect  to  that  of  the  primary 
reflection.  This  indicates  that  for  each  reflection  there  are  more 
sensitive  variations  in  phase  rather  than  in  amplitude  due  to  the 
variation  in  profile  shape. 


1516.5  m/s 
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Fig.  5.6  Geoac'”  '  c  model  for  the  Hatter  as  Abyssal  Plain. 
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RELATIVE  AMPLIT! 


Fig 


.5,7  A  comparison  between  experiment  and  asymptotic  method  using 
the  geoacoustic  model  shewn  in  Fig.  5.6. 
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Fig.  5.9  A  comparison  between  experiment  and  asymptotic  method  using 
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Fig.  5.11  The  reflected  fields  for  the  geoacoustic  model  shown  in  Fig. 
5.6  at  z  +  zq  =  182.7  m. 
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Fig.  5.13  A  comparison  between  the  total  field  and  total  reflected 
field  for  the  c^-l  inear  profile. 


Idien  the  total  reflected  fields  interfere  des  tructively  wi  th  the  direct 

field  at  long  ranges (r  >  2800  m)  as  shown  in  Fig.  5.13,  one  sees  a  rapid 

amplitude  decay  with  horizontal  range  in  the  measurement.  For  the 
2 

1/c  -linear  profile,  the  total  reflected  field  interferes  construc¬ 
tively  with  the  direct  field  at  these  ranges  as  shown  in  Fig.  5.14.  One 
sees  here  that  the  change  in  profile  shape  makes  the  greatest  contribu¬ 
tion  to  the  phase  variations,  and  not  to  the  amplitude,  which  can  result 
in  a  large  difference  in  the  total  field  as  the  reflected  fields  inter¬ 
act  wi  tli  the  direct  field. 


Summary: 

In  both  examples,  the  primary  reflection  is  found  to  be  the  major 
contributor  to  the  total  reflected  field.  In  the  first  case,  the  sound 
velocity  graaient  at  the  bottom  is  large  but  the  value  of  (z  +  zQ )  is 
also  large.  In  the  second  case,  the  sound  velocity  gradient  at  the  bot¬ 
tom  is  small  but  the  value  of  (z  +  zQ)  is  relatively  small.  Due  to 
the  large  water-path  geometric  spreading  in  the  first  case  and  the  small 
souna  velocity  graaient  at  the  bottom  in  the  second  case,  the  multiple 
contributions  become  small.  Also,  the  primary  reflection  coefficient 
amplitude  remains  large  up  to  the  range  of  interest  so  that  the  specular 
reflection  coefficient  amplitude  stays  lower  for  all  ranges  in  both 
cases(cf.  Fig.  3.10  for  example).  Therefore,  in  general ,  a  large  value 
for  (z  ♦  zQ)  enhances  the  contribution  from  the  primary  reflection  at 
1  ong  ranges. 

It  is  also  found  that  the  major  effect  of  the  shape  of  the  profile 
in  the  bottom  layer  is  on  the  phase  variation  of  the  reflected  field 
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which  leads  to  the  characteristic  gross  features  of  the  total  field 
amplitude  as  the  direct  water  field  is  coherently  added  to  the  reflected 
field.  However,  this  effect  would  be  negligible  if  the  water-path  domi¬ 
nated  the  total  geometric  spreading  as  the  source  and  receiver  are 
placed  at  high  elevations  above  the  ocean  bottom. 

Geological  data  near  the  experimental  sites: 

In  order  to  relate  the  physical  features  of  the  geoacoustic  models 
used  in  this  chapter  to  the  actual  geological  properties  of  the  exper¬ 
imental  sites,  seismic  and  core  sample  data  from  the  Deep  Sea  Drilling 
Project(D. S.D.P. )  were  searched.  The  major  objective  was  to  identify 
the  subbottom  reflector  in  our  geoacoustic  models  from  the  D. S.D.P.  data. 

The  seismic  reflection  data  at  60-320  Hz  taken  near  the  D. S.D.P. 
sites  116  and  117,  which  are  the  closest  locations  to  the  experimental 
site,  indicate  the  first  reflector  depth  varying  from  150  to  275  meters 
even  though  their  physical  locations  are  very  close  to  each  other[3]  . 

It  is  therefore  inappropriate  to  extrapolate  the  depth  of  the  first  re¬ 
flector,  namely  the  subbottom  interface,  from  these  data  vrfiich  were  col¬ 
lected  at  least  2.5  degrees  in  latitude  away  from  the  actual  site.  The 
presence  of  a  subbottom  reflector  at  38  meters  below  the  waterbottom 
interface  at  the  experimental  site  was  confirmed  by  3.5  kHz  normal  inci¬ 
dence  records  taken  at  the  site.  Also  yellow  mud,  and  yellow  mud  and 
shells  were  collected  from  19-foot  and  25-foot  core  samples, 
respectively,  at  two  different  locations  within  the  area[l]. 

In  the  case  of  the  Hatteras  Abyssal  Plain  the  depth  of  the  sub¬ 
bottom  interface  agrees  reasonably  well  with  that  of  the  first  acoustic 
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horizon  from  the  drilling  site  387  on  the  western  Bermuda  Rise[4]. 
Horizon  correlates  approximately  with  the  top  of  a  sequence  of  mid¬ 
dle  Eocene  siliceous  turbidites  (lithologic  sub-unit  3A)  which  underlies 
upper  Eocene  radiolarian  mud  (Unit  2).  Other  acoustic  reflectors  were 
also  well  defined  and  their  geological  correlations  coul  d  be  determined 
from  the  core  samples  with  reasonable  confidence.  The  subbottom  depths 
resulting  from  seismic  reflection  and  refraction  profiling  by  Emery  et. 
al  .[5 j  and  Houtz  and  Ewing[6]  near  the  experimental  site  were  130  meters 
ana  215  meters,  respectively,  and  consistent  wi  th  Frisk's  estimate  of 
160  meters.  Therefore,  it  suggests  that  the  subbottom  interface  may 
correspond  to  the  first  acoustic  horizon  AT  observed  near  the  drilling 
site.  The  information  on  geological  properties  of  deep  ocean  sediments 
is  often  not  available;  thus  performing  the  bottom  acoustic  experiments 
near  well  logged  D.S.D.P.  sites  may  increase  confidence  in  geoacoustic 
modeling  and  testing  of  theoretical  results. 
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Chapter  6:  Contributions  and  Related  problems 


In  this  final  chapter  we  shall  summarize  the  results  obtained  in 
Chapters  2  through  5  and  discuss  the  contributions  made  here  to  the 
field  of  ocean  bottom  acoustic  Interaction.  We  shall  also  discuss  other 
related  problems  raised  during  this  research  and  possible  techniques  to 
solve  them. 

6.1  Contributions 

host  of  the  mathematical  discussion  in  this  thesis  was  devoted  to 
the  analytic  evaluation  of  ocean  bottom  reflected  acoustic  fields  given 
in  the  form  of  the  Sommerfeld  Integral .  Within  the  formulation  we  used 
realistic  geoacoustic  models  and  examined  the  sensitivity  of  the  fields 
generated  to  the  geoacoustic  parameters. 

The  asymptotic  method  for  approximating  the  Sommerfeld  integral 
was  derived  in  Chapter  2.  The  total  reflected  fields  were  given  by  a 
series  of  asymptotic  solutions  to  the  fields  associated  with  eigenrays. 
Since  the  reflection  coefficients  as  weight  factors  in  the  integrand 
were  derived  analytically  based  on  realistic  geoacoustic  models  using 
Langer  uniform  asymptotic  approximation,  the  resulting  asymptotic  solu¬ 
tions  were  also  analytic  expressions  for  these  geoacoustic  models. 

The  accuracy  of  the  method  was  confirmed  by  comparison  with  numer¬ 
ical  Hankel  transform  results  except  In  regions  where  critical  incidence 
ano  caustic  formation  take  place.  We  have  also  derived  the  asymptotic 
expressions  based  on  the  Langer  reflection  coefficients  for  the  inter- 


face  wave  propagating  along  the  subbottom  interface  corresponding  to  the 
n-th  multiple  reflection. 

In  Chapter  3,  a  new  canonical  analytic  expression  for  the  sound 
velocity  profile  was  introduced.  We  demonstrated  in  Appendix  D  that  a 

3 

sound  velocity  c  expressed  by  c  -linear  in  depth  provides  a  very  simi¬ 
lar  profile  shape  to  the  one  observed  from  field  measurements.  Using 
the  analytic  expression  for  the  phase  integral  for  this  profile  type,  we 

obtained  complete  analytic  solutions  to  the  reflected  fields  based  on 

2 

this  geoacoustic  profile  shape  as  well  as  for  the  cases  of  the  1/c 
ana  c-1  inear  types. 

We  then  applied  the  above  results  to  examine  several  types  of 
geoacoustic  models  varying  in  layer  thickness,  sound  velocity  gradient, 
profile  shapes  and  source -receiver  geometry  and  demonstrated  complex 
features  of  energy  partitioning  of  eigenrays  with  horizontal  range. 
Because  of  the  accuracy  of  the  plane  wave  reflection  coefficients  we 
aer ived,  our  method  should  be  valid  even  for  geoacoustic  models  with  a 
thin  layer(i.e.  thickness  less  than  a  wavelength)  where  ray  theory  is 
believed  not  to  be  valid. 

The  effect  of  the  sound  velocity  gradient  at  the  bottom  on  the 
geometric  spreading  for  near-grazing  incidence  was  qualitatively  dis¬ 
cussed.  The  conditions  under  which  the  layer-path  of  the  n-th  multiple 
dominates  the  total  geometric  spreading  were  also  derived  for  this 
case.  It  was  proved  that  the  differences  in  profile  shape  were  minimal 
for  a  small  gradient  and/or  a  small  layer  thickness  for  fixed  endpoints 
of  the  sound  velocity  in  the  layer.  However,  in  the  case  of  a  large 
layer  thickness  and  small  gradient,  slight  differences  in  profile  shape 


caused  relatively  large  differences  in  caustic  position  though  the  field 
amplitudes  were  very  similar  including  the  fields  near  the  caustics. 

p 

In  Chapter  4,  we  derived  the  caustic  equations  for  the  1/c  and 
3 

c  -linear  profiles.  We  found  that  the  cusp  of  the  caustic  from  the 
2 

1/c  -linear  profile  always  corresponds  to  a  grazing  angle  of  ir/6  inde¬ 
pendent  of  the  source-receiver  geometry  and  the  sound  velocity  gradi¬ 
ent.  We  also  demonstrated  that  the  leading  term  of  the  asymptotic  ex¬ 
pansion  for  the  fields  near  the  caustic  is  directly  influenced  by  the 
amount  of  curvature  in  the  caustic  curve. 

Two  sets  of  field  measurements  were  discussed  in  Chapter  5.  The 
energy  partitioning  of  the  eigenrays  was  analyzed  and  it  was  found  that 
the  primary  reflection  dominated  the  total  reflected  field  at  most 
ranges  for  both  cases.  We  also  discussed  the  accuracy  obtained  in  mod¬ 
eling  the  ocean  bottom  relative  to  the  geometric  orientation  of  the 
source  ana  receiver  for  given  bottom  structures. 

While  the  density  and  sound  velocity  discontinuities  at  the  bound¬ 
ary  interfaces  determined  the  gross  features  of  the  amplitudes  of  the 
reflected  fields,  sensitive  ray-phase  variation  due  to  the  shape  of  the 
profile  in  the  layer  became  important  as  the  direct  field  was  coherently 
addea  to  the  total  reflected  field.  However,  vrfien  the  source-receiver 
heights  were  large  so  that  the  variation  of  the  water-path  with  respect 
to  that  of  the  layer-path  became  large,  the  total  field  amplitude  became 
insensitive  to  the  profile  shape.  The  discussion  on  this  issue  will  be 
continued  in  the  next  section. 

The  group  of  acoustic  properties  of  the  bottom  which  affect  the 
reflected  fields  given  in  this  thesis  is  clearly  not  complete.  We  have 
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just  started  to  deal  with  mathematical  models  representing  the  basic 
acoustic  properties  of  a  horizontally  stratified  bottom  that  one  may 
find  in  the  real  world.  Without  even  complicating  these  models  there 
are  still  many  problems  remaining.  A  survey  of  such  problems  will  be 
aiscussed  next. 

6  .2  Related  problems 

Additional  aspects  of  acoustic  bottom  interaction  will  be  dis¬ 
cussed  in  this  section.  Some  of  these  are  experimental  subjects,  others 
are  purely  mathematical. 

6.2.1  Suggested  deployment  depths  of  the  receiver  system 
and  an  iterative  inverse  scheme 

We  extend  here  the  discussion  given  in  Chapter  5  and  propose  a 
simple  scheme  to  recover  the  acoustic  properties  of  the  ocean  bottom 
using  a  forward  modeling  technique  such  as  the  one  developed  in  this 
thesis.  We  have  concluded  in  Chaps.  3  and  5  that  in  order  to  model  the 
bottom  accurately,  including  the  profile  shape,  one  should  minimize  the 
sum  of  the  source  and  receiver  heights  to  enhance  the  layer-path  over 
the  water-path . 

This  indicates  that  it  might  be  desirable  to  make  a  preliminary 
estimate  of  the  layer  thickness  by  means  of  standard  vertical  profiling 
(such  as  measuring  the  travel  time  of  a  pulsed  signal  propagating  normal 
to  the  ocean  bottom),  and  then  to  deploy  the  system  of  receivers  accord¬ 
ing  to  the  estimated  layer  thickness  for  oblique  incidence  work.  One 
receiver  should  be  deployed  at  the  depth  where  the  source-receiver  sum 
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does  not  exceed  the  layer  thickness.  The  other  should  probably  be 
located  at  a  higher  elevation  off  the  bottom.  We  realize  from  the  dis¬ 
cussion  in  Sec.  3.2.1  that  the  source-receiver  sum  shoul  d  be  determined 
from  the  sound  velocity  gradient  at  the  bottom,  not  from  the  layer 
thickness.  However,  since  this  quantity  is  not  known  at  this  stage  of 
the  experiment,  estimating  the  lower  end  of  the  sum  of  the  source  and 
receiver  heights  from  the  layer  thickness  is  the  best  one  can  be  expect¬ 
ed  to  do. 

We  have  seen  that  the  fields  measured  at  long  ranges  by  the  top 
receiver  are  dominated  by  the  interaction  between  the  refracting  rays  of 
primary  reflection  and  the  direct  water-path  field  {the  other  types  of 
reflections  are  suppressed  by  their  greater  geometric  spreading).  This 
indicates  that  it  may  be  possible  to  match  the  sum  of  the  direct  and 
primary  reflection  to  the  gross  feature  of  the  real  data  at  long  ranges 
by  using  an  even  simpler  geoacoustic  model  with  a  monotonically  increas¬ 
ing  souno  velocity  profile  (of  the  c-1  inear  type  for  example)  in  the 
lower  halfspace  without  the  subbottom  interface.  As  both  models  (with 
ano  without  the  subbottom)  show  reasonable  agreement,  we  are  obviously 
measuring  most  sensitively  the  sound  velocity  and  sound  velocity  gradi¬ 
ent  at  the  water  bottom  interface.  The  success  Frisk  had  in  obtaining  a 
geoacoustic  model  using  the  parabolic  equation  method[l]  can  be  explain¬ 
ed  by  these  same  arguments,  since  specular  reflection  was  not  explicitly 
included  in  those  calculations  and  the  effect  of  the  subbottom  (they 
modeled  It  as  a  soft  boundary)  was  negligible  at  long  ranges. 

By  matching  the  fields  at  both  receivers  at  short  ranges  with  the 
specular  reflection  included  (but  not  multiples),  one  fixes  the  subbot- 
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tom  interface.  Finally,  by  matching  the  fields  at  long  ranges  at  the 
lower  receiver  with  multiples  now  included,  one  can  "fine-tune"  the  pro¬ 
file  shape  as  well  as  the  sound  velocity  discontinuity  at  both  inter¬ 


faces  . 

Of  course,  the  actual  search  for  geoacoustic  models  which  match 
the  data  is  not  quite  as  strai cfitforward  as  the  procedure  described 
here,  which  is  based  on  purely  mathematical  models.  It  may  also  prove 
impossible  to  suspend  the  source  very  close  to  the  bottom  to  see  layers 
of  snail  thickness. 


6.2.2  Application  of  the  Langer  reflection  coefficient 
to  a  numerical  method 

It  was  noted  that  the  evaluation  of  the  plane  wave  reflection  co¬ 
efficient  via  the  propagator-matrix  method  involves  some  difficulties, 
both  in  selecting  a  step  size  for  the  profile  and  in  computing  time[2]. 
As  shown  in  Fig.  2.2,  the  accuracy  of  the  Langer  reflection  coefficient 
is  quite  good  and  as  long  as  the  geoacoustic  model  is  of  the  type  con¬ 
sidered  in  this  thesis,  the  (total  or  partial)  Langer  reflection  coeffi¬ 
cient  could  be  substituted  in  the  Sommerfeld  integral  to  be  evaluated 
numerically  by  the  Hankel  transform  method[3,4].  This  scheme  should 
also  be  a  good  alternative  to  the  asymptotic  method  introduced  here  when 
it  fails  at  critical  incidence  and  caustics  vhere  more  accurate  fields 
are  desirable. 
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6.2.3  A  partial  expansion  of  the  total  Langer  reflection 
coeffici  ent 

We  expandea  the  total  Langer  reflection  coefficient  before  the 
evaluation  of  the  Sommerfeld  integral.  This  approach  was  suitable  to 
distinguish  one  ray-type  field  contribution  from  the  others.  However, 
it  is  not  necessarily  the  most  efficient  method  for  computing  the  total 
reflected  fields.  Specifically,  the  multiple  contr ibutions  involve 
long,  tedious  summations,  and  we  are  often  only  interested  in  the  total 
contribution  of  the  multiples,  but  not  in  the  separate  multiple  reflec¬ 
tions. 

Therefore,  the  following  expansion  is  suggested  instead: 


vrtiere  the  first  and  second  terms  represent  specular  and  primary  reflec¬ 
tions  respectively,  as  calculated  in  Eq.  (2.27),  but  the  last  term  is  a 
composite  expression  of  the  rest  of  the  multiples. 

The  first  two  terms  may  be  substituted  in  the  Sommerfeld  integral 
to  be  evaluated  asymptotically;  these  provide  ray-type  field  solutions. 
Hcwever,  the  last  term  involves  singularities  for  which  the  normal  mode 
representation  may  prove  ideal.  The  evaluation  of  the  total  contribu¬ 
tion  due  to  the  multiples  would  then  involve  locating  residues  in  the 


complex  p-plane.  This  procedure,  though  it  may  involve  numerical  diffi¬ 
culties,  would  probably  be  more  efficient  than  the  one  used  in  this  the¬ 
sis,  particularly  Wien  a  large  number  of  multiples  contribute  to  the 
fiel  d. 


6.2.4  Limiting  cases  of  interference  headwaves[2] 

In  the  discussion  of  interface  waves  in  Sec.  2.3.3,  we  argued  that 
the  fields  associated  with  those  waves  propagating  along  the  water- bot¬ 
tom  interface  are  due  to  the  saddle  point  contributions  of  the  multi¬ 
ples.  This  is  true  whether  the  sound  velocity  at  the  top  of  the  layer 
is  less  than  or  greater  than  that  of  the  water  column. 

However,  if  the  sound  velocity  profile  in  the  layer  approaches  a 
constant  and  the  layer  medium  becomes  Identically  homogeneous,  the  phys¬ 
ical  characteristics  of  the  interface  waves  become  completely  different 
depending  on  the  sound  velocity  difference  at  the  water-bottom  inter¬ 
face.  Unly  when  Cq  <^(0),  do  they  become  the  classical  head  wave 
associated  with  a  branch  cut  contribution[5 ] . 

It  is  also  of  some  interest  to  show  analytically  that  the  pole- 
-type  singularities  from  the  total  Langer  reflection  coefficient  become 
the  branch  cut  due  to  the  branch  point  at  l/CjiO)  as  c^z)  approach¬ 
es  Cj (0)  for  all  2. 
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Appendix  A:  Accuracy  of  The  Langer  uniform  asymptotic  approximation 


Quantitative  error  analysis  of  the  Langer  uniform  asymptotic 
approximation  is  discussed  in  this  appendix.  We  have  discussed  that 
E(z)G(z)  vanishes  asymptotically  in  Chapter  2.  The  exact  expression  for 
E(z )  is  given  by 


E(z)  J£ 


(A.l) 


where 


n1/6 


(A.  2) 


Thus  E(z)  may  be  further  reduced  to 


EM  .  -  s  qV2  *  a: .  £ 

36  4Q‘  2Q 


(A.  3) 


The  condition  for  the  Langer  solutions  to  be  valid  is 


E(z) 

~T7l 

Id  Q 


«  1 


( A.  4) 


where  the  denominator  is  given  by  the  second  term  in  (2.3). 
Substituting  (A. 3)  in  (A. 4),  we  obtain 


36  4Q*  2QJ 


( A.5 ) 
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Wien  n  is  very  large,  i  .e.  Wien  2  is  away  from  a  turning  point  z^  so  that 


n  s  J  Q(z')  dz'  »  1, 


the  above  condition  may  be  further  reduced  to 

-2r3Q'  Q" 

w  L — tt  “  — <<  1  (A. 6) 

4q  2QJ 

This  is  nothing  but  the  error  criterion  for  the  WKBJ  method[l].  (A. 5) 
and  (A. 6)  are  valid  only  Wien  Q  ^  0  and  n  4  0,  and  Q'  and  Q"  are  not 
large,  i.e.  when  z  is  away  from  z^  in  a  slowly  varying  medium. 

In  order  for  the  Langer  solutions  given  by  Eq.  (2.7)  to  be  valid 
anywhere,  we  must  show  that^never  becomes  zero  even  at  the  turning 
point  although  one  might  expect  such  condition  by  substituting  n  =  0  in 
Eq.  (A. 2).  Let  us  assume  near  the  turning  point  that 

Q2  =  az  as  z  »  0  (A. 7) 


where  zt  =  0  is  assumed  for  convenience  and  a  is  a  constant  gradient. 


Hence  becomes 


(if) 1/6 

which  is  constant  at  the  turning  point,  thus  E(z)  vanishes. 

Therefore,  is  non-zero  anywhere. 

The  Langer  uniform  asymptotic  approximation  is  equivalent  to  the 
WKBJ  approximation  with  an  Airy  function  correction  for  depths  near  the 
turning  point.  However,  when  the  WKBJ  method  is  applied,  it  is 
necessary  to  determine  where  to  switch  from  the  WKBJ  solution  to  the 
Airy  function  solution  in  order  to  obtain  solutions  across  the  turning 
depth.  The  greatest  advantage  of  L.U.A. A.  is  that  it  provides  a  general 
closed  form  solution  valid  anywhere  without  requiring  the  determination 
of  this  boundary  and  continuously  switches  one  solution  type  to  the 
other . 

It  is  also  important  to  note  that  when  Q2(z)  is  linear 
everywhere,  \Jns  constant  and  E(z)  vanishes  entirely,  and  therefore  the 
L.U.A.A.  reduces  to  the  exact  solution  of  (2.1). 
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Appendix  b:  The  Nodi  fled  Method  of  Steepest  Descents  for 

the  Sommer feld  Integral  with  Langer  Reflection  Coefficients 

The  method  of  steepest  descents(M. S.D. )  is  an  expansion  technique 
to  evaluate  asymptotically  an  integral  of  the  form, 

r°° 

I  =  I  F(p )  e1t4i(p)  dp;  u  »  oO  (B.l ) 

-co 

where  F  = 

with  V  =  R  e-2iun‘1 

n  n 

where  is  a  partial  Langer  reflection  coefficient  defined  by  Eq. 

(2.30)  and  u  is  a  layer-phase  defined  by  Eq.  (2.39). 

However,  F  (p )  given  by  Eq.  (2.36)  should  be  written  as 

F  =  F(p,u) 

because  the  modified  reflection  coefficient  Yp  is  a  function  of  p  and 
u.  The  reason  why  we  may  still  apply  the  basic  method  of  steepest  de¬ 
scents  technique  is  that  Vn  behaves  as  if  it  were  a  function  of  p 
alone  when  its  asymptotic  behavior  for  u  >  co  (i.e.  the  dependence  of  F 
on  u  is  very  weak).  Although  we  may  treat  Yn  as  a  function  of  p 
alone,  it  is  still  Important  to  verify  if  V'n  and  V"n  have  the  ap- 
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propriate  behavior  as  a  part  of  the  correction  term  N  in  Eq.  (2.41)  so 
that  N/u  asymptotically  vanishes  as  u  approaches  infinityCl].  For  the 
simplest  example,  let  us  consider  the  specularly  reflected  field  Pq 
and  compare  its  correction  term  N  to  that  of  Lanto  's  problem  studied  by 
Breknovsk  ikh[2] .  Since  both  share  the  same  phase  j6,  we  need  only  con¬ 
sider  the  Rayleigh  plane  wave  reflection  coefficient  given  by  Eq.  (2.31) 
and  the  Unger  reflection  coefficient  in  Eq.  (2.29).  We  have  al¬ 
ready  observed  that  Rq  in  Eq.  (2.29)  may  be  reduced  to  Eq.  (2.31)  when 
the  asymptotic  approximations  of  the  Fock  functions  are  substituted,  but 
this  is  not  the  case  for  their  derivatives.  We  shall  examine  the  asymp¬ 
totic  behavior  of  V'n  and  V"n,  and  justify  the  expansion  necessary 
for  the  method  of  steepest  descents. 

First,  before  we  consider  the  Langer  reflection  coefficient  VQ 
(=  R0),  i.e.  the  plane  wave  reflection  coefficient  for  specular  re¬ 
flection,  we  shall  review  the  leading  order  behavior  of  its  asymptotic 
approximation  which  identically  becomes  a  local  plane  wave  reflection 
coefficient  at  the  ocean  bottom  interface  given  by 


ply0  -  P0Q0 
P1Y0  +  P0Q0 


(B.2) 


whose  first  derivative  with  respect  to  p  is  given  by 


2  PQP  xP 


(B.  3) 
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where  the  subscript  r  indicates  that  this  is  the  Rayleigh  reflection 
coeffi  clent. 


Now  we  differentiate  Vq  with  respect  to  p  and  obtain 


.(2) *  (2)  „ (2)„ (2) 1  v 


which  has  the  leading  order  asymptotic  behavior 


{piYo+poV  Q0 


The  last  term  within  the  parenthesis  is  introduced  in  addition  to  the 
result  from  Eq.  (B.3).  Interestingly ,  when  the  sound  velocity  profile 
in  the  layer  has  the  form  1/ c2 -1  i near,  the  first  and  third  terms  ex¬ 
actly  cancel  out  and  Eq.  f B .5 )  becomes 


2  pqPiP  Q0 

- 2  — 

*piYo  +  poV  V 


(B  .6) 


This  could  be  significantly  different  from  the  value  given  by  Eq. 
(B.3)  because  the  terms  within  the  parenthesis  in  Eq.  (B.3)  tend  to 
cancel  each  other  vhen  cQ  and  c^O)  have  similar  values,  and  the 
aifference  can  become  a  fraction  of  the  order  of  each  term.  In  any  case 
Eq.  (B.6)  Is  still  Independent  of  u  and  behaves  well  in  the  correction 
term  N. 
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ke  shall  next  consider  the  second  derivative  of  Vq, 


V" 

v0 


(A0 


( 2) "g( 2) 


(2)' 


-  W>  -  »i 

b*2)2  ^  ''o21  *  sj,21 


(B.7) 


where  aJ2),  B{QZ) ,  A*2)  ’  and  B^2)'  are  given  by  Eqs.  (2.21-22),  (2.51-52), 

i2)  "  l  2) " 

respectively.  Also  Aq  and  Bq  are  given  by 


aT'--V<tV  wi2>“,IM2l<V 


C0Y0 


♦  (v  v  "  -  v*  )  w,{1),{2)(v  )] 

'Vo  y q  V  w  1 VJ 


(B  .8) 


,( 2)" 


=  -p0L2  S^v0v0'  ♦  SQv’2  ♦  SH  v0v0"]  Wn,»(2)(v0) 


+  frH  c  I  2*1  U|  (  1)  ,  (2)  .  \ 

*-S0  '  S0V0V0  ■*  "  (vg) 


(B  .9) 


Again,  after  substituting  the  asymptotic  approximations  of  the  Fock 
functions  in  Eq.  (B.7),  we  obtain 


V"  - 
v0 


2p0PiP  r#  _  x  Q0  2iplPQ0 

z  [(piyo  poV  zrs - T~ 


(plY0  +  poQor 


Pc  qy  0  Y  0 


2un'(ipgYQ  +  (1  -  i  )p^0q  + 


ip ,  y' 


ro 


)] 


(B.10) 
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where  there  Is  a  term  u  which  leads  to  the  violation  of  the  criterion 
imposed  for  the  asymptotic  expansion  given  by  Eq.  (2.42).  In  an  analo¬ 


gous  manner,  one  can  show  that  asymptotically  Vn '  is  Independent  of  u 
while  Vn"  depends  on  u  for  n  >  0. 

The  direct  cause  of  this  problem  is  due  to  the  mathematical  prop¬ 
erty  of  the  Airy  functions  A  vrfiich  satisfy  the  Airy  equation 


ol^A(x ) 


=  x  A(x) 


ax 


(B.ll ) 


where  x  contains  the  term  dependent  on  u.  As  the  terms  in  Vn  contain¬ 
ing  the  Airy  functions  are  differentiated  twice,  these  quantities  appear 
in  the  asymptotic  form  of  Vn  without  being  canceled  out.  Thus,  unfor¬ 
tunately,  a  direct  application  of  the  M.  S.  D.  fails  at  the  inclusion  of 
the  second  derivative  of  the  Langer  reflection  coefficient. 

However,  we  know  that  Vn  is  a  very  slowly  varying  function  of  p 
ano  it  would  not  make  much  physical  sense  to  expect  a  large  value  from 
the  second  derivative  of  Vn  with  respect  to  p  compared  to  the  first 
derivative.  We  also  know  that  in  practice,  the  major  source  of  limit¬ 
ation  of  the  M.  D.  S.  is  due  to  a  large  value  of  V  vrfiich  is  the  only 
place  where  V"n  could  be  appreciably  large  (ex.  the  Langer  reflection 
coefficient  at  critical  incidence).  All  of  the  above  suggests  that  the 
basic  approach  of  the  M.  S.  D.  is  still  valid.  We,  therefore,  set  V"n 
to  zero,  which  is  equivalent  to  limiting  the  Taylor  expansion  of  F  in 
Eq.  (2.41)  to  somewhat  lower  order  by  rewriting  the  truncated  series  as 


(B.12) 


Ftp) 


F(PS)  +  F'(ps)(p 


1  2>s/2  1 

ps>  +7[v;(-r+77f-» 

2  Y0  ps  Y0 


.  v  1/2  -1„  -2  .  3ps  1 
Vnps  Y0  (2y0  ~T  _  7T)]{p 
Y0  4ps 

We  note  that  F"  In  Eq.  (2.48)  is  the  coefficient  of  (p 
above  expansion. 


ps> 


p$  )2  in  the 
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Appendix  C:  Asymptotic  approximation  of  Fock  functions 


Asymptotic  expansions  of  Fock  functions  are  derived  in  this  appen¬ 
dix  as  follow: 


W(1)*(2)(v)  =  Ai{-v)  ?  i  B1  (-v) 

~  (2v/3irx)1/2e*i(x  '  */4) 
.  v-l/4w-l/2  e±i(x  -  n/4) 


(C.la  ,b) 


( C.2a  ,b) 


W,(1)*(2)(v)  -  »  .-dBi  (-v) 

1  1  ~  d(-v)  d(-v) 


v  *i»/6u(  1),(  2) 


717?  e  . H2/3’^'(x) 


~  -v(2/3irx)1/2e±i  (x  '  3v/4) 


>vl/4  ,-1/2  e*1(x  -  ir/4) 


viler  e 


.  (2  *^/3 


( C.3a  ,b) 


(C.4a,b) 


The  regions  of  validity  of  the  argument  x  are 
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-it  <  arg  x  <  2ir 
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-2ir  <  arg  x  <  * 


for  the  second  kind,  (2) 
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Appendix  0:  Approximation  of  Hamilton's  Profile 

Sound  velocity  profiling  in  ocean  bottom  sediments  from  sonobuoy 
aata  is  given  in  terms  of  the  one-way  travel  time  t  of  a  pulse  signal 
penetrating  through  the  bottom  structures.  It  has  been  found  from 
measurements  that  the  depth  z  may  be  expressed  by  the  third  order 
polynomial  in  t  as[l] 

z  =  At  +  Bt2/2  ♦  Ct3/3.  (D.l) 

By  differentiating  both  sides,  we  obtain  sound  velocity  v  in  t  given  by 
v  a  dz/ot  s  A  +  Bt  +  Ct2.  (D.2) 

Unlike  the  case  where  z  is  a  quadratic  function  of  t[2] ,  it  is  not 
possible  to  eliminate  t  from  the  above  equations  to  obtain  v  as  a 
function  of  z.  Hamilton  does  provide  an  expression  for  sound  velocity 
as  a  function  of  depth  which  is  the  result  of  a  regression  of  points 
generated  by  the  above  equations  for  many  values  of  t  and  is  expressed 
by  third  order  polynomial s[l].  However,  it  is  not  in  a  form  for  which 
the  phase  integral  and  ray  equations  (travel  time,  intensity,  etc.)  can 
be  computed  analytically. 

In  this  appendix  we  shall  Introduce  an  approximate  method  to 
obtain  a  simple  analytic  depth  dependent  expression  for  the  sound 
velocity  from  which  these  analytical  expressions  can  be  obtained 
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easilyL3].  The  results  of  this  simple  velocity  function  are  accurate 
enough  to  model  ocean  bottom  sound  velocity  structure  based  on  the  mea- 
s  uremen  ts . 

Me  first  propose  the  following  quadratic  function  in  z  versus 
souna  veloci ty  V  to  the  third  power: 

V3  =  a  ♦  62  ♦  yz2  (D.3) 

The  above  expression  for  V  has  been  chosen  because  the  shape  of  V 
is  similar  to  that  of  Hamilton's  profile.  We  have  also  assumed  3  un¬ 
known  coefficients,  a,  a  and  y,  to  be  consistent  with  the  3  parameters 
given  in  (D.l)  ana  (D.2)  though  it  will  be  shown  later  that  the  last 
term  on  the  right  hand  side  may  be  dropped  without  sacrificing  much  ac¬ 
curacy.  As  a  result,  the  ray  equations  and  the  phase  integral  can  be 
expresseu  in  simple  analytic  forms  as  given  in  Sec.  3.1. 

After  we  substitute  z  given  by  (D.l)  in  (D.3)  and  retain  terms 
through  t? ,  (D.3)  becomes 

V3  5  a  ♦  8  At  ♦  (8  8/ 2  ♦  y  A2)t2.  (D.4) 

Since  A,  B  ana  C  are  on  the  order  of  unity[l]  and  t  is  typically  less 
than  1  sec,  and  a,  8  and  y  are  also  on  the  order  of  unity  (which  will  be 
verified  shortly),  V  given  by  (D.4)  is  a  good  approximation  to  Eq.  (D.3). 

We  next  cube  both  sides  of  (D.2)  and  compare  the  coefficients  of 
each  order  of  t,  ana  finally  obtain 


8  =  3AB 

Y  =  3(C  ♦  B2/2A). 

From  the  example  given  by  Hamilton  of  terrigenous  sediments  whose 
instantaneous  velocity  in  kn^sec  is  given  by  the  form  of  (D.2)  with  A  * 
1.511,  B  *  1.041  and  C  =  -0.372(Table  III  of  Ref.  1),  we  get  a  =  3.450, 

8  =  8.989,  y  *  1.630.  From  these  nuntoers  it  can  be  easily  shown  that 
the  higher  order  terms  of  t  which  have  been  neglected  in  (D.4)  are  in 
fact  much  smaller  than  the  rest. 

2 

Furthermore,  since  z  is  typically  less  than  1  km,  the  z  -term 
may  be  dropped  in  (D.3)  and  we  finally  obtain 

V3  =  a  +  8Z  (D.5) 

Mien  z  is  .1  km,  V  computed  from  (0.5)  has  less  than  0.2%  error 
compared  to  Hamilton's  formula,  and  even  at  z  =  1  km  where  the  maximum 
error  is  expected,  i  t  has  less  than  1%  error.  Hamilton's  profile  and 
(0.5)  are  plotted  in  Fig.  0.1. 
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